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SYNOPSIS 


The Korringa-Kohn-Rostoker coherent-potential 
approximation { KKR-CPA ) is a first-principles parameter-free 
theory for the electronic structure calculation of purely 
random disordered alloys within the local density 
approximation of the density functional theory. The KKR-CPA 
IS a mean-field theory, in which the random disordered medium 
is replaced by a translationally symmetric effective medium 
with an effective scatterer on each site of the lattice. The 
strength of the effective scatterer is determined by a 
self-consistency condition that the average extra scattering 
(with respect to the medium) from a real atom embedded in the 
effective medium be zero. It has been found that the KKR-CPA 
provides a reliable description of the electronic properties 
of the random alloys. However, the KKR-CPA is a single-site 
approximation, in which the statistical correlations between 
the configuration fluctuations at a given site and the rest 
of the medium are not taken into account. In other words, 
the correlated scattering between two or more sites is 
neglected. Therefore, it is likely to fail whenever a 
single-site effective medium description of the system is not 
adequate, e.g. in presence of strong local environmental 
effects, such as the short-range order and clustering 
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tendencies. The local environmental effects give rise to 
sharp structures in the density of states, which may strongly 
influence the physical properties, like the formation of the 
magnetic moments, the Fermi energy, the density of states at 
the Fermi energy, the electronic specific heat etc. Since 
these local environmental effects involve correlations among 
many sites, they can be investigated accurately and 
convincingly, only within a raultisite or cluster theory. 
This underlines the need for the cluster generalization of 
the KKR-CPA theory. 

The primary aim of this thesis is to present a first- 
principles parameter-free cluster generalization of the 
KKR-CPA theory, within the local density approximation of the 
density functional theory. We have developed this cluster 
generalization by combining the augmented space formalism and 
the conventional KKR Green’s function formulation of the 
disordered alloys. The augmented space formalism was 
originally developed in the tight-binding framework, and it 
provides a self-consistent cluster coherent-potential 
approximation (CCPA). This CCPA preserves the analytic 
properties (herglotz properties) of the averaged Green’s 
function (Chapter 2), and therefore, guarantees a 
single-valued, continuous and non-negative density of states 
at all energies. It also restores the translational 
invariance of the effective medium. These are the striking 
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leaturc'S of our CCPA, which ail other self-consistent cluster 
approaches, with a few exceptions, fail to maintain, and 
therefore, yield unphysical results. The ideas of the 
augmented space formalism are rather general and can be used 
for muffin-tin models as well. We have been able to 
effectively combine the augmented space formalism and the 
conventional KKR method to develop the KKR-CCPA theory. 

In Chapter 2, we present a brief account of the 
augmented space formalism and discuss the unique way in which 
the configuration averaging is handled. The basic philosophy 
behind the CCPA and the procedure leading to it, are briefly 
outlined. One of the achievements of the present work is the 
proof of analyticity of the augmented space CCPA. We prove 
three theorems regarding the analyticity of the configuration 
averaged Green’s function. It is shown that the CCPA yields 
an averaged Green’s function which is herglotz and that, the 
iteration scheme for self-energy always converges to a 
unique, herglotz and bounded value, provided one starts off 
with a bounded and herglotz initial guess. 

In Chapter 3, we present the KKR-CCPA formulation. This 
formulation is valid for purely random systems and does not 
take into account the non-random effects like short-range 
order. Although, in the KKR method there is no off-diagonal 
disorder, the KKR-CCPA, quite predictably, introduces 
off-diagonal corrections as well in the scattering matrices, 
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along with the usual diagonal corrections. A formula for the 
configuration averaged Green’s function within the KKR-CCPA 
is derived from which various electronic properties, such as 
the density of states (total, component, and 1-decomposed ) 
and the charge densities (total and component) can be 
calculated. The component charge densities are essential for 
a full charge self-consistent implementation of the KKR-CCPA 
within the local density approximation of the density 
functional theory. This makes our KKR-CCPA, a 
first-principles parameter-free cluster generalization of the 
KKR-CPA theory of disordered alloys. The problem of 
embedding a cluster of impurities in the KKR-CCPA medium is 
also discussed in detail. The calculation of local densities 
of states on these impurities is helpful in understanding the 
structures in the density of states calculated within the 
KKR-CCPA. 

Chapter 4 deals with the application of the KKR-CCPA 
formulation to a simple model, a muffin-tin alloy in one 
dimension. A one-dimensional muffin-tin alloy retains many 
of the features of its three-dimensional counterpart, while 
greatly reducing the computational effort. Therefore, it 
presents itself as an ideal model for testing the KKR-CCPA 
theory. We have calculated the density of states for this 
model and found that the KKR-CCPA produces sharp structures 
in the density of states in contrast to somewhat smooth 
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representation in the KKR-CPA. We have been able to 
correlate the sharp structures in the KKR-CCPA density of 
states with those in the local density of states on the 
impurities embedded in the KKR-CPA medium. Therefore, we 
conclude that the structures in the KKR-CCPA density of 
states are due to the correlated scattering from two or more 
sites . 

The KKR-CCPA formulation presented in Chapter 3 is valid 
for purely random alloys. However, in reality, almost all 
the disordered alloys exhibit a certain degree of short-range 
order. In Chapter 5, we present a formal extension of the 
KKR-CCPA theory to include the effects of short-range order. 
This theory reduces exactly to the KKR-CCPA formulation for 
random alloys (Chapter 3), in the limit of no short-range 
order. In this formulation, the configuration averaged 
Green’s function in presence of short-range order is expanded 
about that of a random system. This expansion contains an 
infinite number of terms. To make the theory computationally 
feasible, we used an approximation, in which we neglected the 
higher order terms in the expansion. This approximation was 
earlier suggested by Gray and Kaplan [Phys. Rev. B 24^ 1872 
(1981)], in the context of a tight-binding model. We have 
applied this formulation to calculate the density of states 
for the one-dimensional model presented in Chapter 4. It is 
found that, the density of states remains positive throughout 
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the whole energy x-ange, for a reasonably wide range of 
short-range order parameter. It becomes negative at some 
points, only in the strong short-range order regime. This is 
because we have neglected the higher order terms in the 
expansion of the Green’s function, which may be important in 
this regime. Our results are consistent with the results of 
Gray and Kaplan, who used a similar formulation for a 
tight-binding model. 

In Chapter 6, we summarize the main conclusions and 
discuss the achievements, shortcomings, and limitations of 
the present work. Also, we suggest the scope of future work 


in this direction. 
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CHAPTER 1 
INTRODUCTION 


The disordered alloys are an indispensable part of most 
of the present day technological and scientific developments 
taking place around us. Their technological potential can be 
fully exploited if we understand the diverse physical and 
chemical properties of these materials. The electrical, 
magnetic, transport, thermal, optical and superconducting 
properties, the phase formation and the stability of the 
disordered alloys are few examples of such properties. Due 
to their technological importance, the disordered alloys 
have, since long, been an important research topic in 
experimental physics. The theoretical understanding of these 
alloys offers a great intellectual challenge. Most of the 
current theoretical studies on disordered alloys are being 
devoted to the simplest kind of disordered alloys, namely, 
the substitutionally disordered alloys. In such alloys, 
there exists a lattice structure, but each lattice point can 
be occupied by either of the constituent atoms with certain 
probabilities, which depend upon the concentrations of the 
constituents as well as on the local environment of the site. 
In this thesis, we will focus only on the substitutionally 
disordered alloys. 
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The disordered alloys lack the translational symmetry or 
the periodicity, and therefore, Bloch’s theorem is not 
applicable to these systems. With the violation of the 
Bloch’s theorem the entire edifice of elegant theories built 
around it crumbles. This difficulty has been eliminated with 
the development of sophisticated scattering theories and the 
use of Green’s functions. The difficulty arising due to the 
loss of periodicity is overcome by introducing mean-field 
theories. In a mean-field approach, the disordered medium is 
replaced by a translat ionally symmetric effective medium with 
an effective atom on each site of the lattice. Since the 
effective medium is periodic, we can use the existing tools 
for a periodic solid to calculate the properties of 
disordered alloys. 

To understand, how to define the effective medium, let 
us consider the problem of interpretation and comparison of 
properties calculated for the effective medium with the 
experimentally determined results. In an experiment, one 
measures the properties of different samples of the same 
material. In a disordered alloy, the spatial arrangement of 
atoms, in general, varies from sample to sample. Each such 
arrangement is called a configuration or an ensemble. It is 
thus obvious, that different samples of an alloy correspond 
to different configurations. When we talk about the 
characteristics of the samples we measure, we mean an average 
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trend among the different samples, and we are not interested 

in the variations between the samples. Therefore, the 

experimental characteristics of a material is an average 

characteristics of different samples or different 

configurations. For this reason, it is the configuration 

averaged quantity, determined theoretically, which should be 

compared with the experimentally determined characteristics. 

Therefore, the effective medium should be such that the 

effective medium properties are the same as the configuration 

averaged properties of the alloy. This is the philosophy 

behind all the mean - field approximations. A host of 

mean-field approximations have been proposed for disordered 

alloys, such as the virtual crystal approximation^, the 

2 

average t-matrix approximation , the coherent potential 

3-5 

approximation (CPA) etc. 

The current studies on the substitut ionally disordered 

6-8 

alloys are mostly being done within the CPA. The CPA is a 

mean-field theory for purely random alloys, and is a 
single-site approximation. In the CPA, the real disordered 
system is replaced by a translationally symmetric effective 
medium with an effective scatterer on each site of the 
lattice. The strength of the effective scatterer is 
determined through a self-consistency condition, which 
requires that the average extra scattering with respect to 
the medium, from a single real atom embedded in the effective 
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medium be zero. It is generally agreed that the CPA is the 

0 

best single-site approximation. It preserves the analytic 

9 10 

properties of the Green’s function, ’ and has been applied 

very successfully to calculate the electronic density of 

states^ ^ and transport properties'^ of a wide variety of 

disordered alloys. When the CPA is used for a Hamiltonian 

with non-overlapping muffin-tin potentials, the theory is 

called the Korringa-Kohn-Rostoker c ohe rent -potent ial 

7 8 12 

approximation ’ ’ (KKR-CPA). The KKR-CPA theory provides a 

reliable, first-principles parameter-free description of the 

electronic properties of disordered alloys within the local 

13 

density approximation of the density functional 

theory. ’ In the CPA, like all other single-site 

approximations, the statistical correlations between the 
configuration fluctuations at a given site and the rest of 
the medium are neglected. Therefore, the CPA is likely to 
fail whenever a single-site effective medium description 
becomes inadequate, e.g. in presence of strong local 
environmental effects, such as short range order and 
clustering tendencies. The fluctuations in the local 

environment of a site are responsible for sharp structures in 
the density of states, which may strongly influence some 
physical properties, such as the formation of magnetic 
moments, the Fermi energy, the density of states at the Fermi 
energy, the electronic specific heat etc. An illustrative 
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example is the magnetic behavior of the CuNi alloy. It is 

well-known that Ni atoms will carry a magnetic moment only if 

surrounded by a sufficiently large number of other Ni 
1 G 17 

atoms. ’ Also, the experimental value of the critical 

concentration of Ni below which the alloy no longer remains 
1 8"”20 

magnetic ( ^ 47 %) does not match with the theoretical 

1 8 

predictions based on the CPA. In addition, many of the 

disordered alloys are known to exhibit short-range order, 

e.g. CuPd (Ref. 21) and CuPt (Ref. 22) systems. It has been 

found that, the density of states of CuPd alloys calculated 

23 24 

within the KKR-CPA ’ does not compare favorably with 

25 , 26 

experimental results. The disagreement between the 

experimental and calculated results may be due to the 

single-site nature of the KKR-CPA theory, in which local 

27 

environmental effects are not taken into account. These 

local environmental effects, which involve correlations among 

many sites, can be investigated accurately and convincingly, 

only within the context of a cluster or multisite theory. 

The need to go beyond the CPA was realized almost 

immediately after the introduction of the CPA, and since 

then, several attempts have been reported towards such a 

generalization. In the beginning, all attempts towards 

cluster generalization of the CPA were made in connection 

28-50 

with the tight-binding Hamiltonians, because they are 

easy to handle computationally. These approaches can be 
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broadly classified into two categories: self-consistent 

V, 28-46 14, . * * u 47-50 

approaches and non-self-consistent approaches. 

Most of the earlier attempts were to develop a 

self-consistent cluster generalization of the CPA. The 

29 

n-site CPA or CPA(n) of Nickel and Krumhansl , is a 

diagrammatic approach, in which the scattering by all the 

pairs of sites is taken into account. Similar approaches 

30-32 

have also been reported by several other workers. 

51 

Nickel and Butler discovered through numerical computations 
^that these approaches fail to preserve the analytic 
properties of the configuration averaged Green’s function, 
and therefore, lack proper physical behavior, i.e. they give 
negative, discontinuous, and multi-valued density of states 

The molecular coherent-potential 
(MCPA), the travelling cluster 


in some energy regions 
28 


approximation 

45 

approximation (TCA), and the cluster coherent-potential 

46 

approximation (CCPA) developed in the augmented space 

formal ism^^~®^ are the only approaches, which have been 

proved to be analytic, while preserving the conservation laws 
, ^ , , 56-58 

and fundamental sum rules. In these approaches, no 

attempt is made to take into account the scattering by all 

pairs (or triplets). Instead, the statistical correlations 

are assumed to have a finite range, i.e. upto the sites 

within some chosen cluster, and scattering by the sites 

54,55, 58 

inside the cluster is treated exactly. 
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In the MCPA the disordered medium is considered to be 
a collection of disconnected clusters (or molecules) chosen 
so that the entire lattice can be generated by the 
translation of the points in a cluster. The cluster Green’s 
function becomes a matrix and the scalar CPA self-consistency 


condition is generalized to a matrix self-consistency 

condition. The effective medium produced by the MCPA is not 

translationally invariant. In addition, the MCPA is very 

difficult to implement computationally, becoming impractical 

even for nearest-neighbor clusters in three-dimensional 

36-43 

systems. Some authors tried to circumvent this 

difficulty by imposing further approximations on top of the 

MCPA. The idea was to neglect the off-diagonal matrix 

elements of the self-energy completely, thereby making it 

site-diagonal. The self-consistency is imposed only on 

certain diagonal elements of the Green’s functions. Then the 

theory depends only on a single parameter, the site-diagonal 

43 

element of the self-energy. It is clear that this is an ad. 

hoc prescription, but it simplifies the problem to a great 

extent, because, the matrix self-consistency condition of 

MCPA can now be cast into a single scalar self-consistent 

38 

equation. The self-consistent central-site approximation 

was first proposed by Butler and later was used by several 

39-41 . . ... 


other workers. 


In this approximation, 


self-consistency is imposed between the diagonal element of 
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the averaged Green’s function at the center of the cluster, 

and a diagonal element of the Green’s function of the 

effective medium. This is quite unphysical, because, once 

the self-energy is made diagonal, the Green’s function at the 

center of the cluster is almost independent of the external 
4 3 

medium, and therefore, it can not be consistent with that 

of the effective medium. Therefore, it is not surprising 

that, the central-site approximation shows non-analytic 
51 

behavior. This difficulty is removed in self-consistent 

43 

boundary-site approximation, in which the self-consistency 
is required between a diagonal element of the averaged 

Green’s function at a site on the boundary of the cluster, 

and a diagonal element of the Green’s function of the 
effective medium. The boundary-site approximation is 
certainly an improvement on the central-site approximation, 
because a site at the boundary of the cluster is in intimate 
contact with the effective medium, and therefore, the Green’s 
function at a boundary-site will depend on the effective 
medium, and so making it consistent with that of the 
effective medium is reasonable. The boundary-site 

approximation yields non-negative density of states, but its 
analytic ity has not been proved. Moreover, it is an ctd hoc 

prescription and has not been applied to three dimensional 


systems . 
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4 5 57 

The TCA ’ preserves the translational invariance of 
the effective medium, but is a very difficult theory. It is 
derived through a diagrammatic approach. In this case, the 
chosen cluster travels throughout the lattice, but the 
correlated scattering at a given instant is restricted to the 
sites within the cluster. When the cluster-size is increased 
beyond two sites, the diagrammatic analysis becomes 
unmanageable. Also, it is very difficult to implement 
computationally. To our knowledge, it has not been applied 
to any realistic system yet. 

It is now clear that the self-consistent approaches 
which we discussed in the above, are either too difficult to 
implement computationally, or they lack proper physical 
behavior. Therefore, some attempts were made towards non- 
self-consistent generalization as well. The embedded cluster 

47 

method of Gonis and Garland, was the first attempt in this 

direction. In this work, a cluster consisting of a central 

site and its shell of nearest neighbors is embedded in an 

effective medium, which is determined within the CPA. The 

density of states is then obtained from the site-diagonal 

element of the Green’s function at the center of the cluster. 

The effective medium, however, is not determined 

self-consistently with respect to the cluster. Gonis and 

50 

Freeman have proposed a similar technique to include the 
effects of short-range order. 
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The tight-binding Hamiltonians possess severe limitations. 
Various physical quantities that enter into the calculation, 
e.g. the band centers and the band-widths, are determined 
from a parameterization procedure. Moreover, when the 
band-widths of the different constituents of the alloy are 
different, a tight-binding description requires the 
introduction of off-diagonal disorder as well. The 
parameters that describe off-diagonal disorder can be 
determined only by means of over-simplified interpolation 
schemes. Such limitations are removed when the system is 


described by a muffin-tin Hamiltonian. 


This realization 


encouraged a few attempts towards the cluster generalization 

59-62 

of the CPA for muffin-tin Hamiltonians, and this is the 


main subject of this thesis. The work of Gonis ei al 


59,60 


until recently was the only cluster approach in the 

conventional KKR formulation for disordered alloys. This 
59,60 

work is basically an extension of the work on 

. . _ . , _ , , 47 


tight-binding Hamiltonians by Gonis and Garland. 


In this 


work, a cluster consisting of a central site and its shell of 
nearest neighbors is embedded in an effective medium 
determined within the KKR-CPA method. The density of states 
is then obtained from the site-diagonal element of the 
Green’s function at the center of the cluster. In this case 


also, the effective medium is not determined 
self-consistently with respect to the cluster. It was then 
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expected that the effects of lack of cluster self-consistency 

, 27 

would be small. However, Stefanou al have shown that, 
at least for CuPd systems, the effects of cluster 
self-consistency might be significant. 

The primary aim of this thesis is to present a first- 
principles parameter-free self-consistent cluster 

generalization of the KKR-CPA formulation, within the local 

.. 13 — 15 

density approximation of the density functional theory. 

Before embarking upon the cluster generalization, let us 

summarize the requirements for a theory to give physically 

acceptable results. These are as follows: 

(i) The theory must yield analytic, physically 

meaningful results i.e. it must give single-valued and 

non-negative density of states and momentum spectral 

densities. Analyticity is necessary for the fundamental sum 

rules to be satisfied. The analyticity can be assured by 

requiring that the configuration averaged Green’s function 

63 

and the self-energy be herglotz functions. A complex 

operator function F{z) is said to be herglotz if it satisfies 
the following set of properties: 

(a) F(z) is an analytic function of z in the entire cut 
plane Im z 0. Its spectrum is bounded and F(z) behaves 
like 1/z as z * ± oo. 
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(b) F(z*) = f’*’(z), (1.1) 

where F (z) is the hermitian adjoint of F(z), and 

(c) Im F(z) 5 0, for Im z>0. (1,2) 

The imaginary part of an operator function is defined 
as 

Im F(z) = CF(z) ~ F^{z))/2i (1.3) 

(ii) The theory should preserve the symmetries of the 
underlying lattice i.e. the effective medium should be 
translationally invariant. This is essential for the 
calculation of meaningful momentum spectral densities. 

(iii) It should be applicable to general multi-component 
alloys . 

(iv) It should incorporate multisite correlation effects 
for random as well as non-random alloys. 

(v) It should be computationally feasible. 

(vi) It should become exact in the weak scattering and 
dilute alloy limits. 

(vii) It should incorporate various kinds of disorder, 
such as disorder in the diagonal as well as off-diagonal 
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elements of the alloy Hamiltonian in a tight-binding 
description. 

(viii) It should produce density of states which becomes 
exact as the cluster size increases. 

(ix) It should allow the calculation of one- and two- 
particle properties in essentially the same formalism. 


52-54 

We have used the augmented space formalism to 
develop the cluster generalization of the KKR-CPA theory. 
The augmented space formalism is a method of configuration 
averaging a function of a set of random variables. It was 


originally developed in the tight-binding framework. 


It 


46 , 54,58 

provides a self-consistent approximation, the CCPA , 

28 45 

which is a much simpler theory, compared to MCPA and TCA, 
and is a systematic way of going beyond the CPA. In this 
method, the disordered medium in real Hilbert space is 
transformed to a larger Hilbert space, called the augmented 
space, and the CCPA introduced therein. In the CCPA, the 
effective medium is determined by the self-consistency 
requirement that the extra scattering (with respect to the 
medium) from a real cluster embedded in the medium, averaged 
over all the possible configurations of the cluster, be zero. 
The difference between the augmented space CCPA and other 
cluster approaches lies in the sequence of averaging and 
approximation procedures. In the augmented space CCPA, the 
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averaging precedes the approximation, while in the other 

methods, averaging is done after the approximation. The 

exact nature and the philosophy of the augmented space CCPA 

will be discussed in detail in Chapter 2 (section 2.4). This 

CCPA gives a translat ionally invariant effective medium, and 

, 58 

yields a herglotz approximated averaged Green s function. 

Also, the self-consistent equations for self-energy always 

converge to a unique, herglotz, and bounded value, provided 

58 

we start off with a bounded and herglotz initial guess. 

This guarantees a unique and non-negative density of states 

at all energies. The augmented space CCPA has been used 

extensively for electronic structure calculations in the 

46,54 

tight-binding framework. 

The ideas of the augmented space formalism are rather 

general, and can be incorporated within the conventional 

64 

first-principles KKR methods as well. In the KKR method, 

unlike the tight-binding method, there is no off-diagonal 

disorder. Therefore, the application of the augmented space 

formalism to the KKR method is mathematically much simpler 

than its application to the tight-binding method. However, 

computationally, the KKR method would be much more demanding. 

64 , 65 

The works of Hooker jee and Rajput et at were the first 

attempts towards incorporation of the CCPA within the KKR 

method. In these works, an expression for the averaged 

6 12 

density of states using Lloyd’s formula ’ was developed. 
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This formulation has several limitations. Firstly, it has 

been observed that, use of Lloyd’s formula to calculate 

certain electronic properties of alloys, such as component 

densities of states, 1-decomposed densities of states, and 

Bloch spectral density functions is not appropriate, even in 
12 , 66 

the CPA. These quantities become negative in certain 

energy ranges, if calculated using the Lloyd’s formula. It 

is now generally acknowledged that these quantities should be 

1 2 

calculated by the Green’s function method only. Secondly, 

the Lloyd’s formula does not allow the calculation of charge 
densities. The component charge densities are essential for 
full charge self-consistent calculations within the local 
density approximation of the density functional theory. The 
charge densities can be calculated only by the Green’s 
function formalism. 

In the present thesis, we present a self-consistent 

6 7 

field KKR-CCPA formulation for disordered alloys, which has 
been developed by effectively combining the augmented space 
formalism and the conventional KKR Green’s function 
formulation. The KKR-CCPA yields a translationally invariant 
effective medium. However, unlike the KKR-CPA medium, in the 
KKR-CCPA, the effective scatterers are no longer 
site-diagonal. The KKR-CCPA, as expected, introduces 

site-off-diagonal corrections as well, along with the usual 
site-diagonal corrections in the scattering matrices. Since 
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the theory attempts to calculate the configuration averaged 

Green’s function, it allows the calculation of the charge 

densities, and therefore, in this theory, one can achieve 

67 

charge self-consistency within the local density 

approximation of the density functional theory. Thus, the 

KKR-CCPA formulation presented in this thesis, is a 

first-principles parameter-free theory of disordered alloys. 

The most important feature of this theory is its 

computational feasibility. The KKR-CCPA self-consistent 

equations are presented in such a form, that one has to deal 

2 

with the matrices of rank (i+1) only, where I is the maximum 

number of angular momentum states used in evaluating the 

phase-shifts. For most practical purposes, I = 2 is quite 

satisfactory for d-band metals, which means that, one has to 

deal with (9x9) matrices. Therefore, solving KKR-CCPA 

equations on modern computers will not be difficult at all. 

The work presented in Ref. 67 deals with purely random binary 

alloys. However, the extension of this theory to general 

6 8 

multi-component alloys is straightforward. 

We have applied the KKR-CCPA formulation to a simple 

model, a muffin-tin alloy in one dimension. A 

one-dimensional muffin-tin alloy retains many of the features 

of its three-dimensional counterpart, while greatly reducing 

the computational effort, because the involved k-space 

69 

integrations are evaluated analytically. 


Thus, it is an 
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ideal model to test our theory. We have calculated the 
density of states for this model using the KKR-CCPA. We 

observed a few sharp structures in the KKR-CCPA density of 
states in contrast to a nearly smooth KKR-CPA density of 
states. It is found that the sharp structures in the 
KKR-CCPA density of states are due to correlated scattering 
from clusters of atoms. 

The KKR-CCPA formulation (Chapter 3) is valid for purely 

random alloys. However, in reality, almost all the 

disordered alloys exhibit a certain degree of short-range 

order. In Chapter 5, we present a formal extension of the 

KKR-CCPA theory to include the effects of short-range order. 

This theory reduces exactly to the KKR-CCPA formulation for 

random alloys (Chapter 3), in the limit of no short-range 

order. In this formulation, the configuration averaged 

Green’s function in presence of short-range order is expanded 

about that of a purely random system. This expansion 

contains an infinite number of terms. To make the theory 

computationally feasible, we have used an approximation, in 

which we keep only the first correction term in the 

expansion. This approximation was earlier suggested by Gray 
70 

and Kaplan, in the context of a tight-binding model. This 
approximation is in addition to the CCPA. We have applied 
this formulation to calculate the density of states for a 
one-dimensional muffin-tin alloy. It is found that, the 
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density of states remains non-negative throughout the whole 
energy range, for a fairly wide range of short-range order 
parameter. The density of states becomes negative at some 
energy points, only when the short-range order is very 
strong. The unphysical behavior in the strong short-range 
order regime, is probably due to the breakdown of the 
approximation, in which the higher order corrections are 
neglected . 

The outline of the thesis is as follows. In Chapter 2, 
we give a brief account of the augmented space formalism. In 
section 2.2 we derive the augmented space theorem, which 
gives the formula for the configuration average of a function 
of a set of random variables. In section 2.3, we discuss 
some of the aspects of the disordered alloys, and show how 
the augmented space theorem can be a useful tool in 
determining the physical properties of these systems. Due to 
the large dimensionality of the augmented space for these 
systems, it is not possible to use the augmented space 
theorem directly for computational purposes, and one has to 
resort to some approximations. We have adopted a mean-field 
approximation, the CCPA, which is presented in section 2.4. 
The CCPA yields an approximated averaged Green’s function, 
which is herglotz. In section 2.5, we prove three theorems 
in regards to the analyticity of the approximated averaged 
Green’s function obtained within the CCPA. We prove that, 
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the self-consistent iteration scheme for the self-energy 
always converges to a unique, bounded, and herglotz value, 
provided we start off with a bounded and herglotz initial 
guess. This guarantees that the averaged Green’s function 
must also be herglotz. In section 2.6, while concluding the 
discussion, we outline the merits and shortcomings of the 
augmented space formalism and the CCPA. 

Chapter 3 presents the main subject of the work 
presented in this thesis, i.e. the KKR-CCPA. The KKR-CCPA is 
the product of an effective combination of conventional KKR 
Green’s function method and the augmented space CCPA. In 
section 3.2, we present a brief review of the KKR Green’s 
function method. Various equations in this section are 
presented in a form that lend themselves quite naturally for 
cluster generalization, presented in section 3.3. In section 
3.4, we derive the formula for the configuration averaged 
Green’s function within the KKR-CCPA. Also, in this section, 
we show how various electronic properties of the alloy, such 
as the density of states and the charge densities, can be 
calculated. In section 3.5, the problem of embedding a 
cluster of real impurities in the KKR-CCPA effective medium 
and the determination of local properties are discussed. In 
section 3.6 we make the concluding remarks about contents of 
Chapter 3. 

Chapter 4 deals with the application of the KKR-CCPA 
theory to calculate the electronic properties of a 
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one-dimensional muffin-tin alloy. We describe the 
one-dimensional muffin-tin alloy model in section 4.2. We 
have implemented the KKR-CCPA on this model for a two-atom 
cluster. The KKR-CCPA equations for one- and two-atom 
clusters are presented in section 4.3. The KKR-CCPA equation 
for a one-atom cluster quite predictably reduces exactly to 
the familiar KKR-CPA equation. In section 4.4, we present 
the results of our calculations and their inferences. In 
section 4.5 we give the conclusions. 

In Chapter 5, we present the formal extension of the 
KKR-CCPA theory for the systems having short-range order. In 
section 5.2, we discuss some mathematical preliminaries, 
which we will need in subsequent discussions. In section 
5.3, we derive an expression for the configuration averaged 
Green’s function for the systems having short-range order. 
When the short-range order is limitted to nearest neighbor 
atoms only, it is called a Markovian type short-range order. 
We describe a Markov chain in section 5.4. In section 5.5, 
we present our results. While concluding in section 5.6, we 
enlist the shortcomings and the limitations of this 
formalism. 


Finally , 

in Chapter 6 , 

we present 

our conclusions. 

In 

this chapter. 

the achievements, shortcomings. 

and 

the 

limitations of 

the present 

work are 

discussed . 

We 

also 


discuss the scope of future work on the subject. 



CHAPTER 2 

THE AUGMENTED SPACE FORMALISM 


2.1 INTRODUCTION 

In Chapter 1, we have highlighted the need of 

configuration averaging in the study of the properties of 

52-55 

disordered alloys. The augmented space formalism is an 

straightforward and convenient method to determine the 
configuration averaged properties of disordered systems. 
This approach centers around a new technique for 
configuration averaging a function of independent random 
variables. The essential feature of this approach is the 
unique way, in which the configuration averages are handled. 
The random problem in the real Hilbert space is transformed 
into a large new Hilbert space, called the augmented space. 
The augmented space is a direct product of the real Hilbert 
space and the configuration space, which describes the 
configurations of the system. Thus, the augmented space 
contains all the informations about the real disordered 
system as well as its possible configurations. Therefore, 
the configuration averaging is not a further extra process, 
as in other effective medium approaches. 
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A particular configuration of the disordered system can 
be described by a sequence of random occupation variables 
{n^} in a particular order. A certain probability density 
is associated with a particular configuration. Let 

us consider a physical observable of that particular 

configuration. By definition, the configuration average of 
AGn^}) is given by 

<A«n.»> = I AC{n.}) Wn.O 0 dOi • <2-l) 

We begin with the simplification that various ® 

statistically independent. This condition corresponds to a 
purely random system, and is not valid for disordered alloys 
with short-range order. A formal extension of the theory to 
include the effect of short-range order in disordered alloys 
will be presented in Chapter 5. In the present case, the 
composite probability density 5p({n^}) can be decomposed into 
a product of probability densities of individual random 
variables, i.e. 

iPC^n.O = n P.(n. ). 
i 


( 2 . 2 ) 
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Now Eq. (2.1) can be written as 


<Aan.D> 


AC{n.O n P. (n. ) dn. . 

X V 1 1 1 

1 


(2.3) 


In section 2.2, we will show that this integral can be 
evaluated exactly by the augmented space formalism. 

The outline of this chapter is as follows. In section 
2.2, we derive the augmented space theorem, which in effect 
is a formula for the configuration average of a function of a 
set of random variables. This is applicable to any random 
disordered system. In section 2.3, we briefly describe the 
usefulness of the augmented space theorem in calculating the 
configuration averaged properties of disordered alloys. In 
section 2.4, we present the cluster coherent-potential 
approximation ( CCPA ) for a tight-binding Hamiltonian. The 
CCPA is a generalized mean-field approximation, but the 
nature and philosophy of this approximation is distinctly 
different from other mean-field approximations based on the 
embedding of a cluster in the effective medium. The exact 
nature of the CCPA is also discussed in this section. In 
section 2.5, we prove three theorems regarding the analytic 
properties of the augmented space CCPA. We show that the 
configuration averaged Green’s function and the self-energy 
obtained by this approximation, are herglotz. This 
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guarantees a single-valued and non-negative density of 
states. In section 2.6, v?e enlist the achievements and 
shortcomings of the method. 

2.2 THE AUGMENTED SPACE THEOREM 


The basic idea behind the augmented space formalism 

71 

originated from the recursion method of Haydock. The 

recursion method can be summarized in the following few 
lines. A self-adjoint Hamiltonian in a given basis can be 
transformed into a tri-diagonal Hamiltonian in some new 
basis. Then the local density of states, related to the 
Hamiltonian H as 

p_(E) = - Im < i| (E I - H y^\i > (2.4) 


can be expressed as a continued fraction expansion as 


P. (E) 



n 


Im 



(2.5) 


with a down the diagonal and down the off-diagonal 

n n 

positions of the new tri-diagonal Hamiltonian. The local 
density of states has the following properties. 
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P.(E) > 0, 

and (2.6) 

p.(E) dE is finite. 

1 

Now we note that the probability density P.(n. ) shares 

1 1 

with the density of states, the following properties, 

P^(n^) Si 0, 

and (2.7) 

1 P.(n.) dn. = 1 (finite). 

J 1 1 

This suggests that the probability density could 

be expressed as the imaginary part of a herglotz function and 

. 72 . 

can be written as a continued fraction expansion, i.e. 


P. (n. ) 
1 1 



( 2 . 8 ) 


By the inverse process from Eq. (2.5) to Eq. (2.4), 

P (n. ) can be represented as a matrix element of the 
11 ^ 

resolvent of a self-adjoint operator M (Ref. 52), i.e. 


P. (n. ) 
1 1 


n 


Im < V 


O' 


CnJ - >. 


(2.9) 
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Here, M is a tri-diagonal matrix in the configuration 
space (p . spanned by j= 0, 1, ...(m-1), where m is the 

X J| 

number of possible values n. can take* Thus, cardinality of 

1 

the basis is (m-1). The state > is called the 

ground configuration state and it is defined such that the 
matrix element < i-’q | y is the first coefficient 

in the continued fraction expansion in Eq. (2.8). We may 
mention that the operator may not be unique. However, the 
non-uniqueness o f does not pose any problem in the 

mathematical formulation of the problem. The point, which 
does matter is, whether the continued fraction expansion in 
Eq. (2.8) contains finite number of terms, that is to say, 
whether the operator is of a finite rank. This can be 
assured, if we assume that P.(n.) has finite moments to all 
orders i.e. 

[ P.(n.) n^ dn. is finite for I = 0,1,2, (2.10) 

J 1 1 1 1 

This condition is satisfied by almost all the physically 

valid probability distributions, the only exception being the 

52 

Lorentzian distribution. 
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With the help of Eq. (2.9), we can rewrite Eq. (2.3) as 


<A> = Aan.D 


— Im <t-j;l(n.I - dn. 

n ^O'^i ' I 0^ i 


( 2 . 11 ) 


We can evaluate this integral by integrating out one 

variable at a time. Let us first integrate out the j-th 

variable and call it 5^'. 

J 

= I A«n.«[^- iM <yoK’'/ - 


Let {J'^y be the eigenvectors of with eigenvalues /J 


|fi^> 


(2.13) 


Then, the completeness condition is 




= 1 , 


(2.14) 


^ll^2^ ^ *^12 * 
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Rewriting Eq. (2.12) as 



r 

- ^ 

1 


ACCn.D 

L- Im ^ 

n 

v. 


II 


.J ,,J 


Cn.I - M-’ 




dn^ l^o' 


= <!' 


olj 


<pi 


— y y — T” 

" J J Cn. - ) 

U'^ U'^ J 1 

'-‘i 


= <v 


AGn^n 


— y — 3“ 

" Cn, 


<P' 


We 


now use the identity 


n 


Im 


X + iO - a 


= 6(x - a) 


to get 


f, = <y^l[ AC{n » 5 'Hnj - ><^''’1 


= <‘'ol 1 


And since |n'*> are the eigenvectors 


Of M‘ 


1 

dn^|vQ>. 

(2.15) 

(2.16) 

(2.17 ) 

we immediately 


get 
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^ l"o>- (2.18) 

Note that, AC'In^ >i’^j}-» is the same function of 

as AC{n. ,i?^j)-,n.) is of n. (Ref. 52). If we integrate all 

1 j j 

the variables in the above procedure, after integration, all 

the n. ’s will be replaced by the corresponding M'^’s, i.e. 

3 


<A«n.»> = <f] 


where 


I 1 2 3 

= I'-'o ‘'o '■’o 


0 


•> = n h‘o> 


(2.19) 


( 2 . 20 ) 


The state [ f)’ is the ground configuration state in the 

space i = n , which contains the information about all the 

i i 

possible configurations of the system. ACi^^ }■) is an 

operator functional in the augmented space (P = SIf ® ^ , where 

is the real Hilbert space. AOW^}) is the same function of 

as AGn.}) is of {n.}. Eq. (2.19) is called the 
^ 5^ 54 

Augmented space theorem. ’ A general configuration state 

T 58 

spanning ® is of the form 




V 


m. 




V 





>. 


(2.21) 


This state can be 

points O' = ■{ j ; m . 0)- , 

3 


uniquely described by the set of 
at which there are excitations, and 
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the set of cardinalities C = -{m .}• at those points. For a 
binary distribution (as in binary alloys), the cardinality of 
any excitation is 1 (discussed in detail in the next section) 
and thus the labelling C may be omitted in this case. In 


this notation the ground configuration state is O 



where 0 denotes the null set. We shall also write <f Q f 

^ O'* ' i? 

-eye? _ 

as Q for any operator Q in the augmented space. In this 


notation the augmented space theorem can be written as 


<:Aan.D> 



( 2 . 22 ) 


2.3 APPLICATION TO DISORDERED ALLOYS 

The substitutionally disordered alloys are the simplest 
kind of disordered systems. In such alloys, if we neglect the 
correlation between the sites then, each site can be occupied 
by m different kinds of atoms, for an m-component alloy. 

Let us consider an alloy, described by a random 
tight-binding Hamiltonian, 

H = J^i I i><i I + 2 j iXj I * (2.23) 

ij 

ii^j) 


1 
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The states ■{ | iX s^re Wannier states associated with the 
sites {i}. The Green’s function of the system is defined as 


GCz) = Cz I - H ) . 


(2.24) 


where z is a complex number. We will always work on the 
upper half-plane, Im z > 0, throughout this thesis. Now, Eq. 
(2.24) can be rewritten as 



z z 


00 


E 

n=0 


1_ 

n+1 

z 



(2.25) 


so that, any element G- -Cz) can be written as 

ij 


6. .(z) 



<i 


6 

ij 


z 






3 

z 


+ 


(2.26) 


It follows from Eq. (2.26) that, in the plane Im z ^ 0, 

G. .(z) has no singularity, and therefore, is analytic. The 
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singularities of G, Cz) are thus restricted to the real axis. 
Also, from Eq. (2.26) it follows that, 

6 . . 

G. (z) ~ as z > ±00. (2.27) 

1 Jl z 

We can also show, by simple algebra, that, 

Im 6Cz) = - G*^Cz) G(z) Im z, (2.28) 

which implies that, Im GCz) < 0, for Im z > 0 . 

Since H is hermitian and bounded, it follows from Eqs. 

(2.27) and (2.28) that, 6(z) satisfies the herglotz 
63 

properties defined in Chapter 1. The configuration average 
of Eq. (2.28) gives us 

Im <GCz)> = - GCz)!> Im z. (2.29) 

It is evident from Eq. (2.29) that, the averaged Green’s 
function retains the herglotz properties of the exact Green’s 
function. If we define an effective medium, described by a 
translationally invariant effective Hamiltonian £ (also 
called the self-energy), then the Green’s function for the 
system is 

-1 


J^Cz) = Cz I - E ) 


(2.30) 
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thus giving 


In. JTCz) = - Im (z I - J: } (2.31) 


And since is herglotz, it immediately follows from 

^ 10,58 

Eq . (2.31) that also is a herglotz function of z. 

Conversely, if JX^) is herglotz, then J^(z) is also herglotz. 

For any approximation to be physically acceptable, the 

approximated averaged Green’s function and the self-energy, 

must, therefore, be herglotz. In section 2.5, we will show 

that the CCPA presented in section 2.4, conforms to this 

criteria . 

Let us apply the augmented space theorem to find the 

configuration averaged Green’s function for the system 

described by the random Hamiltonian of Eq. (2.23). For this 

purpose we will have to express H ®^s a function of a set of 

random variables {n.}. We consider a binary alloy A B , for 

1 X y 

which the occupation of a site i can be described by a random 

variable n. defined as 
1 


n . = 
1 


1, when site i is occupied by A 


0, when site i is occupied by B 


(2.32 ) 


The random parameters and in Eq. (2.23) can be 

written in terms of the random variables {n.} as 

1 
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£. =: + (c — £ ^ n « 

1 B ^ A B ^ 1 ’ 


and 


{2.33 ) 


V. = + (V.. + )n. n.+ 

ij BB AA BB AB BA a j 


CV.„ - ) n. + ) n. . 

AB BB ^ 1 BA BB ^ j 


Substituting Eq. (2.33) in Eq . (2.23), we get j-j a 
function of random variables {n^} i.e. H = > for which 
the Green’s function, 


D j ^ (2.34) 

is also a random function. The probability density 
can be written as 


GGn.D = z I - HGn. 


P.{n.) = x4(n. - 1) + y<5(n.), (2.35) 
11 1 1 

which satisfies Eqs . (2.7) and (2.10). Therefore, the 
continued fraction expansion of will converge after a 
finite number of terms. Using Eq . (2.16) we can rewrite Eq. 
(2.35) as 
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P. ( n. ) 
1 1 


n 


Im 


f — ^ — ♦ — 1 

I n. - 1 n. J 

L 1 1 J 


n 


Im 


xy 


( 2 . 36 ) 


n . - X - • 

1 n. - y 


This shows that the operator M for a binary 
distribution is a tri-diagonal matrix of rank 2 and is given 

K 

by 


M 


I (xy) 


(xy) 


1/2 ^ 


1/2 


( 2 . 37 ) 


in a space spanned by [^‘q^ and It is easily seen 

that 


and 


M |^.^> = X 


/ ^1/2 , i. 

+ (xy) 

11 i 1 iv 

M = y 


, .1/2 , i. 

+ (xy) r’o^‘ 


( 2 . 38 ) 


A general state |f3 in the total configuration space ^ 

• 58 

can be written as 


f = n U-’ ^ * with n = 
y n I r/ 

1 


" 0, for i e <y 
^ 1, for i e Cf 


( 2 . 39 ) 
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Now applying the augmented space theorem to Eq . (2.34) 

we get 

<Gan.D> = <f| j^z I - HaM^) lo (2.40) 

where is obtained by replacing {n^} in by 


For a ternary system ABC , Mookerjee and Thakur 

X y z 

have shown that if we define the random variables as 

1 , when site i is occupied by A 

0, when site i is occupied by B (2.41) 

-1, when site i is occupied by C 

then the operator is a tri-diagonal matrix of rank 3 with 
the representation 


"i = 


where 



a = x-z, a = 6-{x-z), 



-6 


/?^ = [ x+z-(x-z)^3 t ^2 ~ ^ y-£5(6-x+z)] 


6 = 


[ 


y(x-z)] /[ x+z-(x-z) 



(2.42) 


(2.43) 


with 



37 


The cardinality of the basis is 2 and the single site 
configuration spaces . are spanned by ^ 1 ^ ’ and j f 

. 1 U X 4Li 

v?ith 1 ^’q^ as the ground configuration state. 

We can determine for any multicomponent alloy if we 

know the concentrations of the different constituents in the 

N 

alloy. The dimensionality of the augmented space is NxM , 

where M is the number of species and N is the total number of 

atoms in the alloy. This is quite large for a macroscopic 
23 

system (N ^ 10 ), and therefore, Eq. (2.40) can not be used 

directly to calculate the configuration averaged Green’s 

function. In principle, Eq. (2.40) can be evaluated 

numerically by using the recursion method, and it has been 

done for a one-dimensional model in the tight-binding 
53 

framework. However, for three-dimensional systems, so far 

no attempt has been made in this direction, because it will 

be computationally extremely demanding. In practice, one 

uses the augmented space theorem with some approximation. In 

the next section, we will discuss the approximation we have 

54 , 58 

used, which we call the CCPA. 

2.4 THE CLUSTER COHERENT -POTENTIAL APPROXIMATION CCCPA3 

The CCPA is a mean-field approximation. In this 
approximation, a finite compact cluster of sites, C is chosen 
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out oi {i) in , and the complement of the cluster is 
replaced by an effective medium. It is immaterial which 
particular set of sites we choose. What is important is the 


relative 

position of these 

sites with 

respect to 

one 

another • 

58 

The configuration 

fluctuation 

is 

allowed 

only 

within 

the cluster, and 

therefore , 

the 

conf igurat ion 


averaging is done only over the possible configurations of 
the cluster. The procedure leading to the CCPA involves two 
steps. In the following, we describe these steps. 

Let Cf be the set of all possible vectors in the 
configuration space, i.e. CT = { j f ^ . In the first step, we 
define a subspace spanned by'{|i,f^>, ieC, f^eff} and 
its complement in ’I', • We partition 'I' into these 
two subspaces, and from Eq. (2.40), by applying the partition 
theorem of matrices, we get 

Gi - Hi - Hj, ] ' = [.I, - h] \ (2.44, 

where , 


H = Hi 4 G^ 

(2.45a) 

Hi = H ^1 , 

(2.45b) 

H2 = ^2 H . 

(2.45c ) 
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and , 


^12 ^ ^2 ' 


] 


-1 


(2.45d) 


(2.45e) 


Here is the projection operator in the subspace , 

and can be written as, 


’.■I I 


i;f ><^i;f , 

<y 


ieC f eCP 
{.y 


and , 




i<eC f eCT 

a 


If we write 


H = Hi + a? , 


then we have , 


(2.46) 


(2.47) 


ae' 


c*c 


’ Y~' £70! _a0!’ tCTt’C"’ 

“ ) H- 


.2 ''2 "12 


o(,a’ 


(2.48) 


We shall now introduce the crucial mean-field 
approximation on 3?, in which the subspace 4'- is replaced by 

Zi 

an effective medium. The effective medium is translationally 
symmetric and is defined by the self-energy J], which is the 
effective Hamiltonian. Due to the absence of randomness in 
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subspace » there is no configuration fluctuation in it, 
and therefore, 


and 


where 



G/D 


.ac< 





(2.49) 


In Eq. (2.49), Yin 9 (nxn) and ( mxn ) matrices 
made out of partitioning the self-energy such that the 
cluster has m sites and n=N-m, N being the total number of 
sites. Substituting Eq. (2.49) in Eq. (2.48) one easily 
gets, 




oa' 



02 eL 


(2.50) 


(C ) 

We note that G^CD “ ^ subscript C inside the 

^ (£ ) 

parenthesis implies that, ^ is the Green’s function of the 
effective medium, calculated on the lattice from which the 
cluster C has been removed. Therefore, we write 


ae 


aa ' 




(t ) 



C? 


(2.51) 
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This completes the first step. The rank of the matrices 

II ^ 

and X is mxM , where M is the number of species in the 
alloy , 

In the next step, we again partition into a subspace 
spanned by {|i,f>, i € C }• , and its complement V'g » 
spanned by ■{ j i , f^ > , c? 0, i e £ }. The partition theorem 
gives us the configuration averaged Green’s function in the 
cluster subspace, i.e. 


[<G>]c = [- I 

where , 

H , 



‘l2 ~ ^1 ^ ^2 ’ 


and 


r = I - hp"' - 


*'2 “ ^2 H '’2 • 


(2.52 ) 


(2.53) 


Here p^ is the projection operator in the subspace \f.> 
and can be written as 


Pj = |0<f| , 

and , 

Pn ~ y I f I • 

^2 L ' a o' 

a 0 


(2.54) 
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We note that, 


= |txi| H [ I |f<,Xf^|] = I < 2 - 55 ) 

a ^ (d cr 0 


From Eqs . (2.47) and (2.51), it follows that Eq. (2.55) 

can be rewritten as, 

,acr’ , ,0Cf ’ , , . . ^ 

^12 " ^a0 ^ ). (2.56) 

We know that, > which is defined in in the subspace 

, contains the real cluster, and hence, is independent of 
the self-energy £ (Ref. 58). Thus we can write. 


h 


12 


= K. 


(2.57) 


We can also write 


h 


1 




«00 
+ « 


(2.58) 


We know that H. is the random Hamiltonian in the cluster 

i 00 

subspace. Therefore, by the augmented space theorem, is 

the configuration averaged Hamiltonian of the cluster, which 
is nothing but the virtual crystal Hamiltonian, i.e. 




(2.59) 



43 


Therefore, Eq . (2.52) can now be written as 

[<G> 1 ^ = [z I - - ae - K r® k'*' ] ^ (2.60) 

This is the configuration averaged Green’s function in 
the cluster subspace, which we get from the augmented space 
theorem. This Green’s function is a (mxm) matrix, if we take 
an m-atom cluster. To get the CCPA equations, this Green’s 
function should be equated with the Green’s function of the 
effective medium described by the self-energy J]. The 
effective Green’s function in the cluster subspace can be 
obtained by the partitioning technique, and we get 

= [-1 - E , ' 

= [z I - El - j . ( 2 . 61 ) 

Here, is the self-energy in the cluster subspace, and 
is a (mxm) matrix. The self-consistency condition requires 
that , 

[ s ]c = He • 

which readily gives us the CCPA equation, 

E = ^ K r® k'’ = 3 CD 


(2.63) 
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In Eq . (2.63), we have suppressed the subscript 1 in J]. 

Since is & symmetric matrix with identical diagonal 

elements, the number of independent equations given by Eq. 
(2.63), is m(m-l)/2 + 1. This number can be further reduced 
by exploiting the symmetry of the cluster. 

2.5 ANALYTIC PROPERTIES OF THE CCPA GREEN'S FUNCTION 

The physical properties of a system are intimately 

related to the analytic properties of the Green’s function of 
74 

the system. A physically meaningful approximation must, 
therefore, ensure that the approximated Green’s function 
preserves the analytic properties of the exact Green’s 
function. These analytic properties, which can be derived 
from Eqs . (2.26) to (2.31), are collectively referred to as 

the herglotz properties, and are stated as: 

( i ) 6Cz) is analytic on the entire complex cut-plane 

In z 7^ 0. Its spectrum is bounded, and Q(.z) behaves as 1/z 
as z ► 00 

(ii) Im 6Cz) - 0» Im z > 0, and 

(iii) g\z^ = GCz*)- 

The first two properties together guarantee a 
single-valued, continuous, and non-negative density of 



45 


states. The third condition is called the reality condition, 

because it is related to the fact that the spectrum of H 

(i.e. the singularities of GCz) ) is entirely on the real 

axis. It reflects the causality principle. 

We have seen in section (2.4) that, it is enough to 

prove the herglotzicity of the self-energy to ensure 

that the averaged Green’s function CGCz)]> be herglotz. 

Let us now restrict z in a compact subspace (S of the 

upper half plane ( Im z > 0), such that, whenever z € , the 

virtual crystal Hamiltonian H K (which is independent 

vc 

of J] ) are bounded, i.e. IIHvcB ^ X and m ^ V. 

Let us first state and prove certain identities, which 
we will need in the subsequent discussions. 

f 

(a) IfA=BCB , then Im A and Im C have the same sign. 
Im A = (A- A^)/2i 


= (l/2i) B(C - c’*')b'’’ 

= B (Im Qb’*', (2.64) 


from which the result follows. 
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(b) If Im A > 0, then Im A ^ < 0. 

Im A'^ = (A'^ - CA"V)/2i = I^A"^ (a’’’ - A) (A 

= - A"^ (Im A) CA“V < 0. 

(c) If Im A ^ X > 0, then |jA ^ || - 1/x. 

j|A~^|| = = j^CA’*’A)"^j^^^ = 

= ||Al|“^ = 1^ (Re A)^ + < 1/x. 

Theorem. 2 . 1 

If Im J] < 0 in then Im JfCX]) < 0. 

Proof 

„00 ^ao" „ -1 , 

Let us call Go ~ ^ . Then y = (z 

« -1 
And since, Im J] < 0 and Im z >0, we have Im y > < 

implies that, 

Im 7' < 0. 

Again, from Eq. (2.51) we get, 

Im X = (IB r) Ei2 



( 2.65 ) 


■1 

( 2 . 66 ) 

I - D- 

, which 

( 2.67 ) 


< 0 . 


( 2 . 68 ) 
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Therefore , H = ^ ^ ImH < 0, since is real. 

Also, we have 


and 




+ x 




<ya ‘ 


(O 0 ) 



» 


(2.69) 


which imply that, 

Im = Im af < 0 

and (2.70) 

Im J”’ < 0 . 


From Eq. (2.63) we get 

Im gCD = I- CK PCD k'*’) = K (Im rCD) K''’ < o. ( 2 . 71 ) 

We note that, Eq. (2.63) is solved iteratively as 


Zii) = (2.72) 

Thus, Eq. (2.71) implies that, if Im £(i) < 0, then 

Im J)(i + 1) < 0, i.e. the sign of the imaginary part of the 
self-energy is preserved in each iteration. 

Now let us consider the CCPA Green’s function as given 
by Eq. (2.60). We have 
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Im 


[<6>] = Im (2 I - S) - K Cim r ) K'^ 


(2.73) 


We have seen that for Im J] < 0, Im < 0 and Im J"" < 0, 

which readily proves that,Im ^ therefore, 


Im < 0> 


(2.74) 


in the plane Im z > 0. This condition guarantees a positive 

58 

density of states. 


Theorem. 2 . 2 

If ||Z|| is bounded, so is ||9CE)||. 

Proof We know that Im (zj ~ - x > 0 forzS^ and Im 

< 0. This implies that S 1/x. Hence, from Eq. (2.51) we 

get, 

11*11 ^ IM IIHI IIEl 

< q^/x, (2.75) 

where m ^ q. Since is bounded above by X, we have 

||H|| < X + lael < X + q^/x, (2.76) 


L 
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and so it is also bounded above. Now from Eq. (2.63) we note 
that 


II5CDII ^ + iKf ItrCDII- (2-77) 

-1 

And Im ^ r > 0 (see proof of proposition 2.1), 

implies that yrcDii S l/r, so that we have 

liKDil - X + Y^r, (2.78) 

and thus is bounded. Thus, the self-energy is always bounded 
in the iteration procedure, if the starting value is bounded. 

Theorem. 2 . 3 

The iterative chain £(n+l) = converges to a 

unique solution, provided -Im ♦C(n) 5: x > 0 and -Im Z( n) 2! y > 
0 at each stage of the iteration and E< 0 ) and K are bounded 
to start with. Here the indices inside the parentheses 
indicate the number of iteration in the chain. 

Proof 

The iteration scheme for the self-energy can be written 
from Eq. (2.63) as 


ZM = 3CE<'n-i)) = * K r*"-!' k’*'. 


(2.79) 
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where , 



(2.80) 


at the m-th iteration. We can rewrite Eq. (2.52) as 




(2.81 ) 


Here we have suppressed the subspace indices in and ^ 
for convenience, and these may be restored wherever 
appropriate . 

For the proof, we shall pick only the even steps in the 
58 

iteration, i.e. for m even, 


and , 


where , 


3C(m-l) = J](m-2) ^{m-2) J] (m-3) 


Xim) = E(m-2) !?(m-l) £ (in-3) 


^(m) = (z I - E(m))"\ 


(2.82) 


(2.83) 


(2.84 ) 


Thus our iteration scheme is as follows: 
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E(0) — + ae(i) = E(0) j?(o) zU-i) — * Ed) 


3e(2) = E(o) !^(i) e’’’(-i) E(2) 3e(3) = E(2) ^( 2 ) E^d) 


E(3) > 3e(4) = E(2) 3?(3) E (1)-^ and so on, 
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where, ][](0) = starting self-energy which we 

assume to be bounded and herglotz. 

From Eq. (2.79) we have, 


A][](m) = J](m) - J](m-1) = K - J^(m-2) 


r -1 


= K r(m-l) r (m-2) - P (m-1) P(m-2)K . (2.85) 


With the help of Eq. (2.80), it can be rewritten as 

AE(m) = K PCm-l) ^ at{m-l) - 3e(m-2) jp(m-2)K’’’ 

= K PCm-l) A3e(m-1) P(m-2)K’’'. (2.86) 


By a similar procedure, from Eqs . (2.82) and (2.83), we 


can get 


A3e(m) = af(m) - af(m-l) 


= E(m-2) !^(m-l) A£;(m-1) J?{m-2) j:’’'(m-3). (2.87) 

Since, Im []5(m)] ^ > X > 0 , we can define, Im [!?(m)] ^ = 

2 

(l/oL ) and hence, Eq. (2.87) can be written as 
m 

A3e(m) = a A2](m-1) « ^ « (2.88) 

m-1 m-1 ^ m-z m-z 

with 

Y = t(ra) r^(m-l)/ot . (2.89) 

m ni 

tN F 

t, r IBH ^ ^ 

i\ (A :112J55JL 
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75 

We shall now polar decompose Y as W J where J = 
^ m m m m 

•f 1/2 

(V V ) and W is a partial isometry from the range of J 
m m m ^ ® 

to the range of and thus ||^^Jj = 1* With this we can 


m* 


write Eq . (2.88) as 


LXim) = , I 1 , AS(in-l) a . W |j . 

m-l I m-1 m-l in-2 iii~2 J ni-2 


(2.90) 


^ “"1 2 

In a similar manner, defining Im [ 1 (m)] = ^ 

m 

> 0, we may write Eq . (2.86) as 


A£( 


m 


) = ^ / , Aae(m-l) /5 . w 1 ^ , 

m-l 1 m-l m— 1 m— 2 m-2 J m— i 


= S , ^ , S' 0 ^ 

m-l m-l m-2 


(2.91 ) 


where , 


and , 


9 with S/ = P(m) kV/? 


m 


/V 

m m 


m 


jtf = w'*’ (3 A3t'(m) P , W ^ . 
m m ' m ' m-l m-l 


(2.92) 


With the help of Eqs. (2.90) and (2.91), we can write 
Eq. (2.92) as 

SH = /3 J ^ ^ 1 J«r 

m I m ^ m m-l m-l m-l m-2 

xf<f,0(WJ^/3w l. (2.93) 

m-3 m-2 m-2 m-2 m— 1 



53 


Now , 


ll‘^m-2 °'n.-l Vlil^ - “m-1 ‘^1-2 ^m-2 V 1 i ^ ’ 


(2.94 ) 


2 "t 2 

because |jW || =1. Here, sup <?>|Q Q ] means the maximum 

^ t- 2 

among all the elements Qj<C'^/||pj| for all states of 

the Hilbert space in which Q is defined. In other words, it 
is the upper bound of the operator (i. We observe that, 

C2 k ^2 - C2K-2 = K r''<- 2 ) 

= K r'^(m-2) [im 

= K [- 4 . -r<- 2 )]] 

= K [- Im K'*' = = -Im 3CE(m-2)) 


= -Im ][]( m-1 ) . 


(2.95) 


Thus, from Eq. (2.94) we have 


<r o ^ 

m-2 m 


-1 \-lll " "m-1 

<P 


(2.96) 
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We now define 'j )> = ci f?)> , and then Eq. (2.96) becomes 

m- 1 


||<? o 1 ^ 1 II - sup 

'•‘^m-2 m-1 m-l" 


<:f 1 (- Im £](m-l )) j^:> 


A' <^'| i) \x> 

m-l 


O-i (- Im )) |:t> 


sup 


<.x\ Im (.'$ (m-l )) Ixy 


Now from Eq. (2.79) and (2.84) we observe that 


(2.97) 


and 


where , 


Im E(m- 1 ) = Im 2J ( m-2 ) 


S‘hni-1) = 


2)(m-2 ) = K (m-2 ) K • 


(2.98) 


(2.99) 


Thus, Eq. (2.97) readily gives us 


ll<f « a . W , < 

U‘^m-2 m-l m-l 


<:f| (- Im 2)(m-2)) |:t:> 


sup 

X 


<5:1 Im (G + 0| C-I- 35(»i-2)) !»:> 


VC 


( 2 . 100 ) 


Now, since Im CG^^) ^ - a > 0 and - Im 3) = -ImE— > 
0, the above expression becomes 
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l!f „ c< W 
"■ m-2 m- 1 m- 1 

immediately get 


2 


b/(a+b) , with a and b > 0. Thus we 



Ct 

m- 


1 


W 

m- 


1 


2 



< 


1 . 


( 2 . 101 ) 


In an exactly similar manner, we can get 

W f = \ < 1. (2.102) 

'• m— 1 mm" 2 

Thus, from Eq. (2.93) we have 

IK„II - ^ ^ ^ 

■* Kll ^ 

From Eq. (2.91) we get 


l|Ai:<")|l i IIV2B- 

«|r 

Since, ^ = - Im r(m+l) (from Eq. (2.95)), and T(m+1) 

mm 

is always bounded, if we start from a bounded E( 0 ) ( see 

proposition 2.2), we have - I) > where 7) is some finite 

number. The same is true for an initial bounded starting 
11*^0^ * from Eqs. (2.103) and (2.104), we finally get, 
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iA5;(»)|| < r," |^J|. 


(2.105) 


which implies that 

||z:(p) -EUI|| ^ ri^ 


P-1 

Kll 


0 as p , q > CC' . 


j=q 


Thus 23(m)’s form a Cauchy sequence and, therefore, will 
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converge to a unique value as m > 00. This condition is 

sometimes called the Fixed point theorem. Note that, the 
only assumption, we made was that, the starting value for J] 
is bounded and herglotz. In practice, one usually starts off 
with the CPA value, which conforms to the above assumptions. 
Thus, our proof is complete. 


2.6 CONCLUSION 

In this chapter, we presented the augmented space 
formalism, and derived the CCPA self-consistent equations 
using this formulation. This CCPA restores the translational 
symmetry of the effective medium. We also showed that, our 
CCPA gives an approximated averaged Green’s function, which 
retains the analytic properties of the exact Green’s function 
and that, the self-consistent iteration scheme for the 
self-energy has a fixed point, i.e. it converges to a unique, 
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bounded, and herglotz self-energy if we start off with a 
bounded and herglotz one. This guarantees a non-negative 
density of states, which is a single-valued function of the 
energy . 

The augmented space formalism presented in this chapter 
can not deal with the systems with short-range order. In 
Chapter 5, we will tackle the problem of incorporating the 
short-range order in this formalism. 



CHAPTER 3 

THE KORRINGA-KOHN-ROSTOKER 
CLUSTER COHERENT-POTENTIAL APPROXIMATION 


3.1 INTRODUCTION 


The Korringa-Kohn-Rostoker (KKR) coherent-potential 

6 - 8,12 , 

approximation (CPA) provides a reliable and realistic 

description of band structure and related electronic 

76-83 

properties of disordered alloys. In this method, the 

one-electron alloy potential is constructed of muffin-tin 

potentials, which are non-overlapping spherically symmetric 

potentials centered on the sites of the lattice. In the 

, 84 

beginning, Lloyd s formula was used in all the electronic 

structure calculations for alloys described by a muffin-tin 

Hamiltonian. The Lloyd’s formula gives an expression for the 

density of states for a system that is made up of an array of 

muffin-tin potentials. This expression was manipulated to 

derive the expression for the averaged density of states for 

alloys. However, the component and 1-decomposed density of 

states, and Bloch spectral density function, when calculated 

using the Lloyd’s formula become negative in some energy 

12,66 

regions. Moreover, one can not calculate charge 

densities in the scheme of Lloyd’s formula. The component 
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charge densities are essential for a fully charge 

8 12 

self-consistent ’ implementation of the KKR-CPA within the 

local density approximation of the density functional 

^ 13-15 

theory . 

This realization led to the development of Green’s 
12 

function formalism for disordered alloys. In this 

formalism, one calculates the configuration averaged Green’s 

function, from which various electronic properties of the 

alloy can be calculated. Within this formalism, one can 

calculate the charge densities, and therefore, it can be made 

fully charge self-consistent, within the local density 

approximation of the density functional theory. For this 

reason, the KKR-CPA is rightly called a first-principles 

parameter-free theory of substitutionally disordered alloys. 

Due to the development of fast k -space integration 
85,86 

schemes, it has been possible to use KKR-CPA to 

calculate the electronic properties of a wide range of 

disordered alloys. Now the KKR-CPA has reached the same 

level of sophistication as the KKR band theory of pure 

8 

metals. A relativistic version of the KKR-CPA has also been 
reported . 

In spite of all the successes it achieved, the KKR-CPA 
remains a single-site approximation. In the KKR-CPA, the 
real disordered medium is replaced by a translationally 

symmetric effective medium with an effective scatterer on 
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each site of the lattice. The effective medium is determined 
by a self-consistency condition that, the average extra 
scattering, with respect to the medium, from a single real 
atom embedded in the effective medium be zero. The KKR-CPA 
might fail whenever a single-site effective medium 
description becomes inadequate, e.g. in presence of strong 
local environmental effects, such as short range order and 
clustering tendencies. These local environmental effects can 
be investigated, accurately and convincingly, only through a 
cluster or multisite approximation. 

In Chapter 2, we presented the augmented space 

fcrn.alisn.^'-” and a self-conaistent cluster 

coherent-potential approximation"^^ * (CCPA) in which one can 

go beyond the CPA in a systematic way. In this method the 

effective medium is determined by a self-consistency 

condition that, average extra scattering, with respect to the 

medium, from cluster of real atoms embedded in the effective 

medium be zero. V?e also noted that, the CCPA restores the 

translational symmetry of the effective medium. The ideas of 

the augmented space formalism are rather general and can be 

incorporated within the conventional KKR method as well. We 

have successfully combined these two methods and have 

67 

developed a self-consistent CCPA in the KKR framework, 
which we call KKR-CCPA. Within this formulation, one can 
calculate the charge densities, and therefore, can achieve 
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charge self-consistency within the local density 

approximation of the density functional theory. Therefore, 

our KKR-CCPA formulation is a first-principles parameter-free 

theory of disordered alloys, in which correlated scattering 

from two or more sites are taken into account. The extension 

68 

of this theory to ternary alloys is straightforward. 

The outline of this chapter is as follows. In section 

3.2 we give a brief review of the KKR Green’s function 

formalism. Various equations in this section are presented 

in a form that lend themselves quite naturally to the cluster 

generalization. In section 3.3, we present the KKR-CCPA 

formulation. Although, there is no site-off-diagonal 

disorder in the KKR framework, the KKR-CCPA, quite 

predictably, introduces site-off-diagonal corrections as well 

in the scattering matrices along with the usual site-diagonal 
fi 7 

corrections. In section 3.4, we derive the expression for 

88 

the configuration averaged Green’s function within the 
KKR-CCPA, from which averaged electronic properties of an 
alloy can be calculated. In section 3.5, the problem of 
embedding a cluster of real impurities in the effective 
medium and the determination of local properties, are 
discussed. Finally in section 3.6, we give our conclusions 
about the implementation of the KKR-CCPA. 
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3.2 THE KKR GREEN’S FUNCTION FORMALISM 

Consider a disordered substitutional alloy A B . The 

X y 

one-electron Hamiltonian for a certain configuration of the 
alloy can be written (in atomic units) as, 


H = - ^ i>.(r - gp = - 7^ + ^ 

i i 


where R.’s are the site 

1 

locations 

- 

and are 

1 1 

the 

potentials centered at 

. The potentials 

are 

random, and can be .) 

A 1 

or (r.) 

Ir x' 

depending upon 

the 

occupancy of the i-th site. We assume 

muffin-tin type and do not overlap with 

that . Ct . ) are 

1 1 

each other: 

of 

1 1 11 M 

- rp 

(3 

.2) 


where r is the muffin-tin radius, and ©(x) is the step 
M 

function, which is 1 when x^O, and vanishes for x<0. 

The Green’s function for this system can be written in 

12 

operator notation as, 

G = G + 6 T G » 

o o o 


(3.3) 
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where is the free electron Green’s function, and T is the 
total scattering operator of the system, and is given by 
Dyson equation: 


T = V + V 6^ T , (3.4) 

O 

where V is the total crystal potential. The operator T 
be written as 


1-1 

ij 


(3.5) 


where s^re called scattering-path-operators 

defined as 


8,12 


and 


are 


Here 


<5 + G 

i j o 


2 

k^i 


t" = C I + ), 

1 O 


(3.6) 


(3.7) 


is the t-matrix that describes the scattering from the 

potential on the i-th site. Due to the spherical symmetry of 

the potentials, the single-site t-matrices are diagonal in 

the angular momentum space, and their on-the-energy-shell 

6,8 

matrix elements are given by 
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LL’ 


- (»)■’ exp[ 


sinf' 


(i) 


6 = t: 6 

LL’ I LL’ 


(3.8) 


where L = is a composite index, Jt = and are 

the phase-shifts for the potential 

1 X 

The free electron Green’s function, when projected onto 
the real space takes the form^’^ 


G^Cr, r’) = -i?£^Y^(r) j^(Kr^) h^(Kr^) Y^(r-’) (3.9) 


where r^(r^) is the smaller (greater) of the variables r and 
r’, Yj^(r) is a real spherical harmonic, and j^(x) and hj^(x) 
are the spherical Bessel and outgoing Hankel functions 
respectively . 

We will derive the expression for GC^ from Eq. (3.3) 

for two different cases; the site-diagonal case, when both r 

« • 

and r’ are inside some i-th cell, and the non-site-diagonal 

4 12 

case, when r and r’ are in different cells. Taking the 

projection of Eq. (3.3) onto the real space, for these two 

different cases, we get 


G (?,?') 


I zic? j iii. 


it c?'p - y 

L j ij / . 


zfc?,) J.c?',). 

Xj X Xt ^ 


LL’ 


L 


(3.10) 
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~t I . 

where r and r lie on i-th and j-th cells respectively. The 

• Ib/ « 

functions Zj^C*" and respectively the regular and 

irregular solutions of the SchrOdinger equation for a single 

12 

muffin-tin potential, i.e. 


V^+ t'. (r . ) - E ] 
11 




= 0 , 


(3.11 ) 


and are normalized such that at r . = r , they join smoothly 

1 M 

to 

and 

= Y (r . ) j, («r. ) (3.13) 

L 1 L 1 i 1 

where 

Cj = = - * cot^l^^ + ix. (3.14) 

In Eq. (3.10), are on-the-energy-shell matrix 

LL 

elements of scattering-path operators. These can be obtained 

. , . . 12 

from Eq. (3.6) by simple algebraic manipulation to give, 


where , 



( A"’) 


ij 

LL' ’ 




A 

^ . . 

IJ 



(3.15) 


(3.16) 
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1 -) 


B, j / ^re the real space version of the well-known KKR 
^ 6 8 

structure constants, ’ and are given by 


Kl- = - 


E 


L" 


l-l> -1*^ 

i ^ h, „(*eR. .) Y_,(R. .) (1-6. .) 

LL’ L ij L ij ij 


(3.17) 


where . =1^. - and , are Gaunt factors, given by 
ij j 1 LL’ 


A 


L" 

LL 


, = J Y ,(X) Y (X) Y (X) dfi . 


(3.18) 


The Fourier transform of Eq. (3.17), 




(3.19) 


are related to the KKR structure functions, S. . ,CK) as , 

LL 


Bll.A = * i* . 


(3.20) 


Eqs. (3.15) and (3.16) will play a central role in our 

T i j 

formulation. We observe that, the randomness in comes 

from C^’s only, while , depend completely on the lattice 
structure. Therefore, in the KKR framework, unlike as in 
case of tight-binding formalism, there is no off-diagonal 


disorder. 
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3.3 THE ICKR CLUSTER CPA FORMULATION 


We now apply the augmented space formalism to the KKR 

Green’s function formalism presented in the preceding 

section. We will determine the configuration averaged 

path-operators using this formalism. The first step in this 

direction is to express the randomness in terms of 

some random variables. From Eq . (3.16), the matrix A can be 

67 

written in an operator notation as, 


A 



L ^ 

ij( ) 


|iXj 


(3.21 ) 


In Eq. (3.21) we have suppressed the angular momentum 
indices. We now describe the randomness in A Uy writing C 
in terms of the random variable n^ defined in Eq. (2.32), as 

C^ = C^ + 6C , (3.22) 


where , 


6C = C^ - C®. 


(3.23) 
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Substituting Eg. (3.22) in Eg. (3.21) we get 

A = C® ^|i><il + 6C ^|iXi| ^ |iXJ| 

^ i ij(i?^j) 

= A({n^)^). (3.24) 

Our aim is to find configuration averaged and by 

the augmented space theorem, Eg. (2.19) it is given by, 


where , 


= < i,f| 



j,f >» 


(3.25) 


AGM'd = 


Y b'-* iixji 

/ 4 immmJt LmU 




f e(S 

a 


6C iXil ® 


(3.26) 


The operator and the vectors j f^ are described in 

Chapter 2 (section 2.2). In Eg. (3.26), AOMi) is an 

operator functional in the augmented space 'I', and CT is the 
set of all possible vectors in the configuration space. Now, 
we partition the augmented space 5* into subspaces and 

in the scheme described in section 2.2 of Chapter 2. 
Subspace is spanned by ■{}i;f^, ieC , f ^ e (£}• , where C is 
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the chosen cluster. Subspace 'P is the complement of 'i' in 

i 

'i . If the cluster C contains m sites, then the rank of 4-'^^ is 
m 

mx2 . The crucial mean ^ field approximation now involves 
replacing the conjugate space by an effective medium. 

Since we need the <i , f | . . . | j , f> element (i, jel^ ) , of 

in Eq . (3.25), we partition AGmS) as follows, 


A = 



(3.27) 


where A^ is in the subspace 4'^ and A 2 i® i’^ subspace . 

By applying partition theorem, the inverse of A in^ subspace 
can be written as, 

[A-h, = 

The four constituent matrices of A are given by, 


ieC f eCT ieC 

a 


® M 


ijeC f e® 

(i?»=j) 


- y 


B 


10 


iX j I ® 


I 


(3.29) 
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A 


2 


keC ' kleC ' 

(k^l) 


gkl 

kleC' 
{k?^l ) 



|kXl|] ® ^ |f^f^ 

f ecr 
a 


(3.30) 


Ai2 = -[ Z Z " E Z Z 


ieC keC' 


ieC keC' 


f e(S: 

c 


(3.31 ) 


and 



kl 

In Eqs. (3.30) and (3.31), b are the site-off-diagonal 

corrections in C, the diagonal corrections already contained 

in C . Note that, A, does not contain b^"^, because 

subspace ' 1 '^^ contains the real atoms and their configurations. 

A„ has got these off-diagonal corrections, because subspace 

'4' contains only the effective atoms. This partitioning of 

the augmented space into subspaces 'i'^^and different 

from the physical partition of the lattice into 

28 

non-overlapping clusters of Tsukada scheme. In the present 
scheme, the subspaces and are interacting, and 

therefore, the operator A. 9 » which connects the subspaces 

X _ X 

, , . 1 . 1 k; 

and , must contain the off-diagonal corrections b , in 

order to include correlated scattering between a site inside 
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the cluster and a site outside the cluster.®^ Also, this 

ensures that, the CCPA restores the translational invariance 

i k 

of the effective medium. Removing b from A in Eq. (3.31) 

JL 2 

will lead to the old Tsukada scheme, which is not the same as 

kl 

our CCPA, In general, b will not be zero in KKR-CCPA, 
However, for a one-atom cluster, these terms vanish, thus 
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For convenience, we 


reducing to the correct KKR-CPA limit. 

, , , kl ,-^kl ^ , 

add D to D (which are also off-diagonal in site indices 
kl 

and define R ^ ^ as , 
ef f 


-,kl ^kl , kl 

Beff = B + b 


(3.32) 


Thus Eqs. (3.30) and (3.31) can now be written as, 




E 


B 


kl 

eff 


kXl 


kleC' 

(to*=l) 




f 

a 


f ecr 

a 


and 


A 


12 



(3.33) 


(3.34 ) 


The elements of A in Eq. (3.28), can be found if we can 
evaluate the triple product (A 22 ^2 ^21^ ‘ Since is the 
matrix in the effective medium with the cluster C removed 
from the lattice, we have 


1 _ -p(C) 

2 “ 'eff ' 


(3.35) 



72 


The superscript C inside the parenthesis indicates that 

the path operator matrix of the effective medium with 

the cluster removed from the lattice. It is not possible to 

calculate this quantity, but we will show that it can be 

6 7 

completely eliminated from the computational procedure. 
From Eqs. (3.33), (3.34) and (3.35) we get. 


A A"' A = y y T 

^12 ^2 ''21 l_^ ‘^eff ' 


(C )kl plj I -W . I- 


ijeC kleC 


I 




f eCT 
a 


1 

ijeC 




f ecr 

Cf 


(3.36) 


where 




T (C )kl -.1 j . ■ ^ 

. for 1.0 e C. 


(3.37 ) 


kleC' 


It is clear from Eq. (3.36) that the matrix [Aj 2 ^£^^21^ 
is diagonal in the configuration space ^ . The off-diagonal 
parts in A thus come from Aj^ only. From Eqs. (3.28), (3.29) 
and (3.37) we get, 



^(C® - I iXi 

ieC 


i jeC 

(i^j) 



9 




f ecr 


& 



9 


(3.38) 
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Now we start the second step of the CCPA, in which is 
partitioned into subspaces and such that is spanned 
I i > i s » and its complement y.' is spanned by 
^ 1 ^ ^ ^ i € C} . In this scheme, we partition A as, 


i e C} 

. In 

this 


f : 

si 

A = 

1 1 ; 

[ • • • 

12 


L ^ 

Sf 


^ 21 

2 •' 


(3.39) 


where, is in the subspace y/ and is in subspace V^o • 

^ 1 Z 2 

The ranks of ^ and •si ^ are m and mx(2™-l) respectively. The 
four constituent matrices of A can be written as 


= 

lO<f| A |f><f|, 

(3.40) 


|fXf| A ( ^ 

f (C^0) 

O' 

f ^ ^ r % 

(3.41) 

^2 - 

il 1 WlJ A [ ^ 

f (c^0) f Acy’^0) 

(3.42) 

= 

■"l2 • 

(3.43) 


These matrices can be written explicitly with the help 
of following identities in the augmented space: 


1/2 


fi>, 


I f> = X I f> + (xy) 


( 3 . 44a) 
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j 

fi> = y 1 

f,> 4 

|f>. 

(3.44b) 

1 

IV = " 1 

* 

(xy) ] f 

ij>’ 

(3.44c) 


^4i> = y 

4> 

X f 

+ (xy) 

If •>. 

( 3.44d) 


' 1 J 

' IJ 

' J 

1 

if, .> = X 


, a/2 

+ (xy) 

|fij,;>, (w ^ 1). 

(3.44e) 


and so on. 

The above identities can be written in a general form as, 

= X I f^ + (xy)^'^^| f^,>, 0’’= O' + i, when i e o 

and (3.45) 

i 1/2 

H I f = y |f^ + (xy) o’= o - i, when i e o. 

Therefore, we can write 

<fl m' 1 = X 6.^ . (3.46) 

With the help of Eq. (3.46), we can write Eqs. (3.40) to 
(3.43) as 

= J\c - |i;f><i;f| - ^ li;fXj;fl, 

ieC ijeC 

(i^j) 


(3.47) 
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'^12 ~ ^ ( 3 . 48 ) 

ieC 

•^21 = ^ ^1 > ( 3 . 49 ) 

ieC 

and 

^2 = [ - ae^b |iXi| - ^ cb"' + |iXj|] ® 

ieC i jeC 

( J 


Y. 


iXi| ® ^ 


O' 0 

ie€ 

O 0 

tyV 0 

(3.50) 

where , 

— 




C 

= X C + 

y C » 

(3.51) 

and 

0) 


6C. 

(3.52) 


Now, inverse of A subspace is given by 

[ = [^1 '■ 

Since the ranks of the matrices ‘^12’ *^21 *^2 

small, it is not difficult to find the matrix elements of 
[A”^]^ in Eq. (3.53). FromEqs. (3.48), (3.49), and (3.50) 
we can write 
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•^12 ^2* -^21 = “ u 

1 

“ “ r (3.54) 

where 

r = ^|i;f><i;fj |j;f^><j;fl. (3.55) 

1 jeC 

From Eq. (3.53) we get 

<i;f| A ^1 j;f> = <i;f| - w <o )~^j j;0- (3.56) 

But, from Eq. (3.25), the left-hand side of Eq. (3.56) 
is the configuration averaged path-operator, i.e. 

<T*''> = <i:(| (*^1 - U r " (3.57) 

We note that, both and J”’ are diagonal in |0> so 
that we can write 

<T^'^> = <il - CO Q (O )‘^| j>, (3.58) 

where 

^ =a| ^ lo 

= Vcc - |iXi| - ^ |iXj| (3.59) 

ieC i jeC 

( J ) 
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and 

Q = <f| r |0 = ^li><i;f.l (3.60) 

i jeC 


Now consider a translationally invariant effective 
medium for which the path-operators are given as 


where 



A 


eff 



i j 

ii^j) 


(3.61 ) 

(3.62) 


We need the path-operators of the effective medium in the 
cluster subspace, which we get by partitioning technique , i . e . 

Tefr [ X - *^"5 |k><k| - 

keC 

®eff ^ l>‘> <1|] 

kleC 
(^ 1 ) 

The self-consistency condition requires that the 
averaged path-operators obtained by the augmented space 
formalism should be equal to the path-operators of the 



-1 

lj>. (3.63) 


effective medium, i.e. 
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<jiJ> ^ -j-ij 
‘eff 


(3.64) 


which implies that , 

Q ^ - JE^h liXil 

ieC 

^^iff * I • (3.65) 

i jsC 
(iJ^J) 

The KKR-CCPA equations follow directly from Eq. (3.65), if 
we compare the corresponding matrix elements on both sides, 



Ceff = C - 6) <i;fi| (3.66) 

= 6i <i;f.| 6J, for ij{i?^j)eC . (3.67) 

1 iS J 


We note that, Eqs . (3.66) and (3.67) are self-consistent 

equations, which involve C , and on the 

eff 'L 

right-hand side as well. These equations can be solved 

iteratively; the starting values for C and b^'^ may be 

ei I 


taken as C and 0 

X j 

evident that, 

( i , j e C ) , and C 


respectively. From Eq. (3.63), 
can be expressed in terms of > 


eff ’ 


it 


IS 


Kf? 


1 . e 
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it’ 


i .1 


= S 


13 


Kif^’ 


eff 




€ c 


(3.68 


In Chapter 4, we will show that, for a two-atom cluster, 
Eq . (3.68) can be cast into a very simple expression. The 
path-operators appearing in Eq. (3.68) can be calculated by 
the usual Ic-space integration over the first Brillouin zone. 


Cl 


1 


= I [C -g -(if)] exp[ -ik. ^ . .] dk 

'eff „ 3 ^ eff “eff ij 

8n ■'BZ 


(3.69) 


where R „„(k) is the Fourier transform of R i.e, 
eff eii 




= R(ic) + 2 ^ ®xp [ 

ieC 


(3.70) 


We also note that, 

= C (3.71) 

and, therefore, we need to solve Eq. (3.67) for j>i only. 
This reduces the number of self-consistent coupled equations 
in Eqs. (3.66) and (3.67) to [ m(m-l )/2] +1 . This number can 
be further reduced by exploiting the symmetry of the cluster. 
Note that, all these equations involve matrices of rank 
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(i+1 ) , where I is the maximum number of angular momentum 

states used in evaluating the phase-shifts. Therefore, with 
the modern computers, solving KKR-CCPA equations for a 5-site 
cluster (for which there are 11 coupled equations) will not 
be difficult at all. 

A cluster consisting of a central site and its shell of 

nearest neighbors should be the ideal choice. But it will be 

computationally demanding and so we propose an alternative, 

though approximate method. A nearest-neighbor cluster has 

the rotational symmetry of the lattice. Thus, one can take a 

two-atom cluster consisting of the central site ( site 0) and 

one of its nearest neighbors (say, site 1). The KKR-CCPA 

equations for a two-atom cluster, are given in Chapter 4. By 

solving these equations, one can obtain C „„ and And 

eff 

, Offi 

then, b , for the rest of the nearest neighbors, can be 
found by rotational symmetry, since, b^"^ will transform like 
. These b^^, can then be used in Eq. (3.70). 

Note that, as in the KKR-CPA, the first-order 
correction to the KKR-CCPA inverse t-matrices is of the order 
of (/ = xy C6C)^. This implies that, the KKR-CCPA, like the 
KKR-CPA, is exact in the limit of vanishing concentration of 
either constituent and also, in the limit when the difference 
between the two atoms is small. 
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3.4 THE CONFIGURATION AVERAGED GREEN’S FUNCTION 


The Green’s functions for site-diagonal (SD) and non- 


site-diagonal (NSD) cases can be written from Eq. (3.10) as 


12 


= Tr[ t“ FV. .?',)] - I 


and 




(3.72) 


(3.73) 


where, Tr stands for trace over the angular momentum index L, 
and 


.r') 


1 1 


and 


,r'.) 

^LL''^ i ’ r 


Z.V.) Z.S (?'.), 

14 1 li 1 




(3.74) 


(3.75) 


Configuration averaging Eqs. (3.72) and (3.73), we get 

<6soC?.?')> = Tra’Wj .?;)> - I <zj;c?.)jj;(?'.3>. (3.76» 


and 


NSD ' ’.II V 


(3.77) 


The second term in Eq. (3.76) can be evaluated directly. 
The joint averages appearing in Eqs. (3.76) and (3.77) are of 
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the form <ff U y and <(X y. These can be found by the 

8 9 

augmented space method, if we can express them in terms of 
the random variables n^ defined in Eq. (2.32). We write, 

(3.78) 

and 

i i BR 

U = U "i Ug n. + U 3 , (3.79) 

where , 

<5U = U'^ - U®. (3.80) 


and 


Uj = 


U2 = - u®®. 


U3 = - u®®. 


(A) Determination of <t“ u^> 


(3.81 ) 


Using Eq. (3.78), we can write 

a“ ub = <t“> u® + <T“n.> 6u 

= T°°f U® + <T“n.> 6U. 


(3.82 ) 
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5 2 

By augmented space theorem, we can write 

where AC{M^}-) is given by Eq. (3.26). Using Eq. (3 
get, 

<T^^n^> = X <i;f| A li;f> + <i;f| A 

But, by augmented space theorem, 

<i; f] A ji'.O = <T^^> = . 

so that, we have, 

<T^^n^> = X + (xy)^^^ <i;fl A 

Substituting Eq. (3.84) in Eq. (3.82) we get, 

<T“ uS = T°°f U + <i;f| A -fu. 

where 

U = X + y U®. 


(3.83) 

44a ) , we 

1 i;f ■>. 


(3.84) 

(3.85) 


(3.86) 
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^-1 

The only unknown term in Eq. (3.85) is <i;f| A ji:f.>. 
It can be found within the KKR-CCPA by following the 
procedure outlined in the preceding section. After the first 
partitioning of the augmented space, we get, 


<i ; f 


~-l .-1 

A 1 i;f ■> = <i;f| A | 


(3.87 ) 


where A is given by Eq. (3.38). After the second 

.-1 

partitioning in the scheme of eq. (3.39), A 

..-1 ^ 

falls in the top right-hand side block of A > i.e. if we 

..-1 

partition A as > 


..-1 


A = 



(3.88) 


..-1 

then <i;fj A I i 5 will fall in block, which by 

partition theorem, is given by 






(3.89) 


where the matrices appearing on the right-hand side are given 
by Eqs. (3.47) to (3.50). With the help of Eq. (3.54), we 
can write Eq. (3.89) as 
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.-1 

<i;f I A I i 


A |i;f.> = - <i;f| - u P u j "^2^ I i i 

III [<'■'! K - “T"] 


j , k f f 

a O'* 




(3.90) 


From Eq. (3.48), we have 


<[ j ; f 1^, ^tk;f ,]> = (i) 6 «5 . 6 

(yi 12' jk <y0 kcv' 


(3.91) 


and thus, Eq. (3.90) can be rewritten as 


<i;f| A |i;f.> 


^ |^<i ; f 1 - CO J7’ Cij j j ; f> 


X <j;f^ j I ^ 


i;V]. 


(3.92) 


From Eq. (3.57), we have 


<i;f| - (^T 


j;f> = T 


eff ’ 


(3.93) 


and from Eqs. (3.66) and (3.67) we have 


" “‘‘h - '=efP ■ '^ip] 


(3.94 ) 
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Substituting Eqs. (3.93) and (3.94) in Eq. (3.92) we get 


<i;f| A 


, -1 V* 

“ ■ 2^ b <. 

jeC 


( 3.95 ) 


Thus, from Eq. (3.85) we finally get 


<t“ uS = T°° J U + - Q (6C)' 


' fiu j - 


Y bj* (6Q‘^ 6U. 


eff 


( 3.96 ) 


jeC 


We note that, 

<T“ US i* <T“XUS. (3.97) 

6,64 

which is true within the single-site KKR-CPA as well. In 

Chapter 4, we will show that, if we take a one-atom cluster, 

Eq. (3.96) exactly reduces to that obtained by the restricted 

averaging method within the KKR-CPA. 

Also, we note that, the expression for the joint average 

within the KKR-CCPA given in Ref. 64 (Eq. (16b) of that 

paper) is not correct, because, in that expression the second 

term of Eq. (3.96) is missing. However, within the 

88 

single-site KKR-CPA, it is exactly the same as Eq. (3.96). 
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(B) Determination of 


Using Eqs. (3.79) and (3.81) we can write 


<T1J ^ ^BB ^ + <niT"'’>U2 + <T'‘'’n XJ 


(3.98 


Let us determine the terms in Eq. (3.98) separately. By 
augmented space theorem we have 


<n.T"'’> = <i;f| I^AGmS)] |j:f>. 


(3.99) 


Using (3.44a), we can rewrite Eq. (3.99) as 


= X Tgjj. + (xy)^^^ A | j;f>. 


(3.100) 


In a similar manner, we can write 


<T '^n.> = X + (xy) <i;f A j;f 


(3.101 ) 


Also, we can write 


<nj‘-’n_.> = <i;fl r ><•> | j ; f > 


(3.102) 
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Using Eq. (3.44a), one easily gets, 


y T-ij . 2 T-ij , .1/2 

<niT n^> = X Tgjj + xCxy) ' <i; 


f| A ij;f •> + 


- ^ 1 

A |J;f> + xy <i;f.| A 


(3.103) 


Putting Eqs. (3.100), (3.101), and (3.103) in Eq. 


(3.98) we get. 


<T” = T” U + A I j;fj> (x Uj U3 ) + 


- ^ -I 

A I j;f> (x Uj + U2) + xy A . 


where 


,, 2 ,,AA 2 ,,BB ,,AB , ,BA 

U=xU +yU +xyU +xyU • 


(3.104) 


(3.105 ) 


The element <i;fj A | j ; f will fall on E block of 

3 1 4 

Eq. (3.88), and therefore, it can be written as 


r 1”^ -1 

<i ; f I A I j ; f^> = -<i;f||j^^-t«>J^“J ■'^12 *^2 ^ 

-I I lb 


k,l f f , 

a cf* 




(3.106) 
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Using Eq. (3.91) we get 


<i;f| A = - ^<i;f I a)j lk;f><0 






(3.107 ) 


Using Eq. (3.94) in Eq. (3.107) one finally gets 


<i;f| A 


J;f > = CC - C) 

/ ' eff eff ^ 



k 


(3.108) 


In a similar manner we can obtain 


^-1 


ij 


<i;f| A =“ ^eK 


y u'l 


eff 

(3.109) 


.-1 


The element <(i;f^ A 


j;f .> falls on E block of Eq. 

j z z 


(3.88), which is obtained by partition theorem, as 


^22 = < - Si S2 


(3.110) 
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Thus, we have, 


<i;f I A 

X ^ 




21 


k,l f f , 

c a ’ 




(3.111 ) 


Using Eqs. (3.67) and (3.91) and a little algebraic 
manipulation we get. 


^-1 


A <i;f j ; f> (C - a 

- V (3.112) 


-1 


From Eq. (3.89) we get 




- y ..-1 b “ 


eff 


(3.113) 


Therefore, all the terms of Eqs. (3.104) are determined 

through Eqs. (3.107), (3.108), (3.112), and (3.113). The 

88 

configuration averaged Green’s function, for site-diagonal 
and non-site-diagonal cases can be obtained from Eqs. (3.96) 
and (3.104) by substituting >^p »*"p for 

and U ^ . 
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The configuration averaged density of states is given by 


<n(E)> = Im 


n 


a 


<Ss/i ""i 


(3.114) 


where the integration is over the central (i-th) Wigner-Seitz 
cell. The configuration averaged charge density is 


<p{r)> 





Im 






•+ 

dr . 


(3.115) 


The site-diagonal form of the configuration averaged 
88 

Green’s function is given by, 


<Gsj,Cr,r')> = Tr 




) + - C) (iQ'Vc?,?’) 


E 


<5FC?.?’) - 


J [x Z*(?) !*(?’) + y Z®(?) (3.116) 

L 


where 


and 


FCr.r’) = X F C^.r’) + y F (I'.r'’), 
(5FCr,r’) = F (r ,r’) - F (r ,r’). 


(3.117) 


(3.118) 
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In calculating the local properties, such as the 
component densities of states and component charge densities, 
one needs the restricted averaged Green’s function. The 
component charge densities are essential for a charge 
self-consistent implementation of the KKR-CCPA, within the 
local density approximation of the density functional theory. 
The restricted averaged Green’s function can be determined by 
embedding a real atom in the effective medium. We will 
discuss the problem of embedding in the next section. 


3.5 EMBEDDING A REAL CLUSTER IN THE EFFECTIVE MEDIUM 


Several authors have discussed the problem of embedding 

47,60,90,91 

a cluster of real atoms in an effective medium. 

, 60 

Gonis et al have discussed how to embed a real cluster in 
an effective medium, determined within the KKR-CPA method. 
The generalization of this technique to an effective medium 

the KKR-CCPA medium is not 
straightforward,”' because, in the KKR-CCPA, the effective 
t-matrices are no longer site-diagonal. We have, for 
convenience, clubbed the off-site-diagonal parts of the 
effective t-matrices ) with and have defined a new 


determined within 

67 


variable B 


ij 

eff 


In the discussion that follows Eq. (3.31), 


we have stressed that, when we embed a cluster of real atoms 
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in the effective medium, the off-diagonal term between a site 

inside the cluster (a real atom) and an effective atom, 

should be so as to include correlated scattering between 

67 

these two sites. Therefore, in the present situation, we 
assume that, the off-diagonal term between two real atoms be 

i j 

B , whereas, that between two sites, when at least one of 

i j 

them is an effective site, be B „„ . Admittedly, this is an 

ef f 

ad hoc prescription and needs further work. In Chapter 4, we 
will discuss the effect of different choices of the 
off-diagonal terms (between real and effective atoms) on the 
density of states. 

In the present prescription, the restricted averaged 
Green’s functions, following Ref. 12, can be written as 


= Tr <t“> 


i=:Ol 


F (r. 


.?p 


1 


(3.119) 


and 




j 


--ft 


Tr 




i=ayj=ft 


r-ftOL, . 

F C*' 




(3.120) 


where a, /? = A or B. The restricted averaged path-operators 

^ . 12 ^ 

can be found by a simple partitioning technique, to give 




=a 


_0i —oo 
^ * eff 


> 


(3.121) 
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and 


<T1J> 

^ 1 , 

jij || 
j=p ^ 'eff 

"eff®'" 



- C )D- t"" 1 

^ efr^ ' ef fj 

“1 -t 




(3.122) 

where 




D' = 

I + - c )1 

_ ^ ' eff ^ eff 

-1 

( 3.123) 


The component density of states for an c-type atom is 
given by 

[d“ T:°, FV„ .?„) d?„] . (3.124) 


(a ) ■+ 

The component charge density p (T ) on an 0!-type atom 


is given by 


92 


P 


“(?) 


1 m X 

Tr Im 

n 



oo 

eff 


dE], 


(3.125) 


One can construct the potentials from p \i' ) using 
the local density approximation of the density functional 
theory, and hence, can achieve charge self-consistency. 
The methodology is as follows. One starts off with the 
potentials of the pure a-type solid ( a = A or B), and 
then solves the KKR-CCPA equations given by Eqs . (3.66) and 
(3.67). From Eq. (3.125), one calculates the component 
charge densities, which are used to calculate new potentials 
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• These new potentials are again used to solve the 
KKR-CCPA equations. This procedure is continued till the 
charge densities in the successive iterations become 
consistent. Therefore, this method is independent of initial 
parameters . 

Note that, if we average Eq. (3.124) over the 

configurations of the n-th site, it would not be identical to 

88 

Eq. (3.114). This is expected, because, Eq. (3.114) is 
derived by averaging over all the configurations of the 
cluster. For a one-atom cluster, both these equations, 
though derived through quite different techniques, are 
exactly identical. We may stress that, while calculating the 
averaged properties of an alloy within the KKR-CCPA, Eqs . 
(3.114) and (3.115) should be used rather than averaging Eqs. 
(3.124) and (3.125). However, Eqs. (3.124) and (3.125) are, 
nevertheless important, when one needs the local properties, 
such as the component density of states and component charge 
densities . 

If we embed two impurities of type ot and ft on the sites 
i and j respectively, the local density of states on the 
respective sites is given by 


and 


n- (E) = 
"2 '"> = 


Tr Im 

Tr Im 


^i=a,j=ft 



(3.126) 


n 


(3.127) 
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Hence, we need only to find the path operator matrices 

][ i J J 

• It is straightforward to 
obtain these matrices by the partitioning technique. The 
real space is divided into two sub— spaces; one for the 
impurity sites and the other for the rest effective medium. 
Then by simple matrix algebra one easily gets 


<t“>. , . ... 

1=C.X, J|=fj 



OO _ _ij 

eff" eff 




+ 


and 




i=£:x, j=ft 



(3.128) 




- C 




(3.129) 


Eqs. (3.124) to (3.129) can be used to calculate local 
density of states and charge densities for a single impurity 
or two impurities of a foreign kind as well. 


3.6 CONCLUSION 

In this chapter we showed that, one can go beyond the 
CPA in the conventional KKR Green’s function formalism by the 
application of augmented space formalism to the KKR Green’s 
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function. We found that, the KKR-CCPA introduces site-off- 
diagonal corrections as well in the scattering matrices, 
along with the usual site-diagonal corrections, which we do 
not get in a single-site approximation like the KKR-CPA. The 
KKR-CCPA restores the translational symmetry of the effective 
medium. The KKR-CCPA formulation has been put on the same 
rigorous footing as the KKR-CPA, in which one can calculate 
the density of states, the component density of states, and 
the charge densities, self-consistently within the local 
density approximation of the density functional theory. 
Therefore, our KKR-CCPA formulation is a first-principles 
parameter-free theory of disordered alloys, which includes 
correlated scattering from clusters of atoms. In this 
formulation, we have neglected the short-range order. A 
formal extension of the present formulation to include the 
effects of short-range order will be presented in Chapter 5 . 



CHAPTER 4 

APPLICATION OF THE KORRINGA-KOHN-ROSTOKER CLUSTER 
COHERENT-POTENTIAL APPROXIMATION TO A MODEL SYSTEM 


4.1 INTRODUCTION 

The implementation of the Korringa-Kohn-Rostoker cluster 

coherent-potential approximation (KKR-CCPA) to realistic 

systems is involved, and is computationally demanding. 

Therefore, it should be first tested on simple models before 

a full-fledged realistic implementation is attempted. This 

formulation has already been applied to an over-simplified s 

65 

phase-shift model. This model rests on the od. hoc 

assumption that, the shape of the density of states curve is 

semi-circular. In this chapter, our aim is to apply this 

formulation to a one-dimensional muffin-tin alloy, for which 

67 

we do not have to make such assumption. Another advantage 

of choosing this model is that, for this model, the involved 

69 

k-space integrations can be obtained analytically, which 
greatly reduces the computational effort. 

In section 4.2, we describe the one-dimensional muffin- 
tin alloy model. We have calculated the density of states 
for this model using the KKR-CCPA formulation for one-atom 
and two-atom clusters. In section 4.3, we outline the 
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methodology of implementing KKR-CCPA for one-atom and two- 

atom clusters respectively. It is found that, for a one-atom 

cluster, the KKR-CCPA equations reduce exactly to the 

67 

familiar KKR-CPA equation, which is the correct limit. 
Also, we show that, the expression for the configuration 
averaged Green’s function for a one-atom cluster reduces 
exactly to that, what one would get by the restricted 
averaging method in the KKR-CPA. In section 4.4, we present 
the results of our calculations of density of states within 
the KKR-CCPA. We observe that, the KKR-CCPA density of 
states is somewhat structured in comparison to the KKR-CPA 
density of states. We show that, the structures in the 
KKR-CCPA density of states are close to the structures in the 
local density of states on the impurities embedded in the 


KKR-CPA effective medium. This suggests that, the structures 
in the KKR-CCPA density of states are due to correlated 


scattering from clusters of atoms. Also, we present the 

averaged density of states on the central site of a cluster 

of atoms embedded in the KKR-CPA medium, for clusters of 2, 

3, and 5 atoms respectively. These calculations are similar 

.60 ^ 


to those proposed by Gonis dl . 


We find that, this 


averaged density of states reproduces some of the structures 
of the KKR-CCPA density of states, which become prominent as 
we increase the cluster-size. In section 4.5, we give our 


conclusions . 
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4.2 THE ONE-DIMENSIONAL MODEL 


We consider a one-dimensional alloy of non-overlapping 

69 

muffin-tin potentials described by the Hamiltonian 



n 



(4.1) 


In close analogy with the muffin-tin model in three 

dimensions we assume that w (x) is symmetric in x and 

n 

vanishes for jxli r (the muffin-tin radius), r S a/2 (a is 

the lattice constant). Due to this symmetry, various 

quantities, such as the t-matrices and wave functions, can be 

69 

expanded in terms of symmetric and anti-symmetric functions 
defined as: 

Y (x) = l/y2, 
o 

and (4.2) 

Yj^(x) = l/y2 sign(x). 


These functions are analogous to the spherical harmonics 

in three dimensions. The one-dimensional analog of the 

69 

spherical Bessel and Hankel functions are defined as. 


j^(z) = cos ( z-ln:/2 ) , 1=0, 1, 

and (4.3) 

h^{z) = exp[ iz-i7r/2 ] , 1=0, 1. 
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Almost all of the familiar identities involving 
spherical Bessel functions and spherical harmonics have 

one-dimensional analogs. 

The on-the-energy-shell matrix elements of the single 
scatterer t-matrix are given by 

~ ^ exp[i0^] sinB^ 6^^, (4.4) 

where « = Ve and 6^ are the phase-shifts. Unlike the case in 
three dimensions, in this case, we have only two components 
of I, i.e. 0 and 1. The site-diagonal form of the Green’s 
function of this system in real space can be written as, 

6Cx,x' ) = J Z^(x) T°°, (X' ) - J z^ix) (x' ), (4.5) 

U ’ L 

where Zj(x) and (x) are respectively the regular and 

irregular solutions of the SchrOdinger equation for an 
isolated muffin-tin potential, and are normalized such that, 
at r = jxj = * they join smoothly to 

Z^(x) = Y^(x) h^(»tr)j, (4.6) 

and 

J^(x) = Y^(x) j^(Jer), (4.7) 

where 

Cj = l/t; ■ 


(4.8) 
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In Eq . (4.5), are on-the-energy-shell matrix 
elements of the path-operator, and for the KKR-CCPA effective 
medium, it can be written as. 


n/a 


T 


1 J a 

eff “ 2n 


J dk C 


-n/a 




*“eff~ ^eff^^^ ik(i-j)a] , 


where 


01 


= B(k) + b exp(ika). 


(4.9) 


(4.10) 


B(k) 

dimension 


are the KKR structure function 
whose elements are given 


matrix 


in 


one 


~ I ®^PC“ik(i-j )a), (4.11) 

i^j 

where B^j, are the real space versions of Bi 2 ^,(k), and are 

,“69 
given by, 


B 


ij 


expCik |i-j| a) 


lit 


-i sign(i-j) 


[ i sign(i-j) 


(4.12) 


The indexing of the sites is as follows: the central 
site is taken as 0-the site, the numbering on the right side 
is in terms of positive integers, and on the left it is in 
terms of negative integers. Substituting Eq. (4.12) in Eq. 
(4.11) we get , 
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B(k) 


i sin(?fa) sin( ka ) 

cos ( ka )-cos(^^a ) cos ( ka ) -cos ( Wa ) 

-sin( ka ) i sin(xa) 

cos(ka)-cos{«a) cos(ka)-cos{Ka) ^ 

<4.13) 


The inverse of the effective t-matrices, C may not be 

ef f 

diagonal in angular momentum space, when calculated within 
67 

Therefore, we consider a general case, when 


the KKR-CCPA, 


are not diagonal in /-space. 


Eq . (4.9) can be 


transformed into a contour integration on a complex z-plane 

67 

by putting z = exp(Lka): 


-1 j 
eff 




dz z 


j-i~l 


g(2) 


4 3 2 

X 2 +6 z +0( z +/3 z +y 


(4.14) 


where C is the unit circle centered at z=0 and g(z) is a 
matrix in i -space, whose elements are given by 


1- 

2z 

cos(Xa) ) 

[i- 

J£b^^zl + 

2z 

sin(xa) 

1 - 

,2z 

cos()ea) ) 



2z 

sin( Xa ) 


(4.15) 

2 2 

§ 01 ( 2 )= ^01 cos(Jfa))+ z - 1 , 


2 2 

giO<z)= ^10 "-^ 10 ^^ cos(»a))- z + 1. 


and 
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Also, we have 


^'"'<*’00 ‘'ll - ‘’oi ‘-10 >’ 


6 = 7) - K - b^^) - 2 X cos(Ka), 


(X = X + ^ + -2 7) cos(«ea) , 


f3 = 2 cos(xa) [2 - r - ^ + 2sin(sea) 


<4.16) 


and 


[ 2i - * . „ . * 

^ tSo Si - Si SoJ - (So ^ Si>- * (Sr So5 


T) = n«Cb +b)-ECCb +Cb-Cb -Cb^, 
' ^00 11^ '^ 00 11 11 00 01 10 10 or 


In Eqs. (4.15) and (4.16), C = C and b 

et I 


wOl 

b ; the 


subscripts are the angular momentum indices. 

Now let us consider a system, where the one-dimensional 
muffin-tin potentials are of the form 


iJ (x) = 0(r - r). 

n u M 


(4.17) 


These are potential wells, and the phase-shifts for a 

93 

single potential can be written as, 
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<y,= - Tt r + tan 

L M 


^ -I 


where 




k^= y E + 


The wave functions Z^(x) inside the muffin-tin 

normalized such that they match with Eq. (4.6) at r = 
93 

are given by 


cos(ier 6 ) 

2 ix) = T—z r; r~ cos(i«r + 6 ) Y (x), 

o K sino cos(k,r^,) o o 

O 1 M 


sin{Xr + 0.) 

ZiO** = sin(kr V ""i'"’- 

1 In 


The density of states can be found as, 


1 _oo r- 

n(E) = - -:: 7 - Im Tr T 2. > ' L t ’ ’ 


n 


where 


.a/ 2 


IV "" .. 


Z, (x) Z, / U) dx 


■a/2 


ll‘ 


2E sin^0 


I 


r cos {Hr^+e^-ln/2) sin(2k^r^' 

___ Ir^ + — 

cos (k.r -tfT/2) 1 

1 M 


a 

2 


H-(-l) 


sin(*<a+20j^ ) - sin{2itr^+2d 


( 


2« 


H 


(4.18) 


cell, 


(4,19) 


(4.20) 


(4.21) 


In) ^ 

—y 


(4.22) 
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The of f-di atonal elements of F vanish due to opposite 

6 7 

parity of 1=0 and 1=1 wave functions. 

The configuration averaged density of states can be 

obtained by using Eq . (3.114), and the local density of 

states can be obtained using Eq. (3.124). However, if we 

assume that, the off-diagonal term between a real atom and an 

67 

effect.iv<* atom is of the form, 

CB'i. - B^'’), 0 < AJ < 1, (4.23) 

new eff 


then Eq . (3.124) can not be used to calculate the component 
density of states. The case /J = l , represents the one which we 
have adopted in deriving Eq. (3.124). In the present case, 
the restricted averaged path-operator appearing in Eq. 
(3.119) can be written as 


= (X + P W ‘ S) ‘|t:“ * P 


i=0( 




(4.24) 


where 


X = 1 cc“- 


P = T°°f 

W = I <1-,.). 

s = Cc“- ,i-„. 


(4.25) 


and 



XU/ 


4.3 KKR-CCPA FOR ONE- AND TWO-ATOM CLUSTERS 


Let UK coriKider a one-atom cluster consisting of the 


0-th 

site* 

The 

conf i guration 

spaniied by 

1 


and 1 f 
‘ 0 

for t 

his case, 

we have 


If 

(C® - 


lo><ol® 

i f ><f 1 

* 0 0 * 

where 






space (J)^ for this site is 
I O'. 

= Rewriting Eq. (3.50) 


+ 6C |0><0lejf Xf |M°|f Xf I 

I ' ' O O' ' o o' 

(4.26) 




00 

1 


jki^O 


r^O ,3 “¥“(0)^1^ 

^ef f ' eff 


B 


ko 

eff 


(4.27) 


With the help of Eq. (3.44b), Eq, (4.26) reduces to 


= (C - 3E°°) |0X0| ® |f„Xf^|. (■'•28) 

where 

C = y + X C® (4.29) 

It is evident from Eq. (4.28) that, <i;f^|.Jf 2 = 0, 

for i J. Therefore, for a one— atom cluster, the KKR— CCPA 
equations given by Eq. (3.66) and (3.67) now become 


and 


■“eff 


= C - 0) (C 





(4.30) 


0 . 


(4.31) 
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Thus, predictably, there is no off-diagonal corrections 
in this case. We note that, for a translationally symmetric 
effect. iv<‘ medium, the site-diagonal path operator can be 
olitained by a partitioning technique to give, 

pOj ..^(o)jk _ko 1-1 
^eff ' eff ^effj 

jk^O 

Eq. (4.32) readily gives us, 



3e 


00 

1 



CT°° ) 

efr 



(4.33) 


Putting Eq. (4.33) in Eq. (4.30) we get. 


"ef f 


0 


^ 1 - 

^eff * ^^eff^ 




(4.34) 


which readily reduces to the 

. 8,12 

KKR-CPA equation, which is 


more 


familiar form of the 


'eff 


= C + (C 


CC® 


- =efP- 


(4.35) 


The site-diagonal form of the configuration averaged 
Green’s function, for this case, can be rewritten from Eq 
(3.116) as 
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<G,„(r,r ■)> = Tr[T"“^ j Rr .r’) * ^ ) (iC)*' 

6F(r ,I ’)|j - ^ [x zjc"-) J^Cf ’) + y Z®(r) J®Cr o] • (4.36) 

L 

By simple algebraic manipulation one can show that, Eq. 
(4.36) can be rewritten as, 

<G 5 ^(r.r')> = Tr X D'' T°“j rV.rO + y D® 7°°^ FV.r’)] - 

I [x Z^cr) J*cr ’) 4 y Z®Cr) • (4.37) 

L 

A B 

where D and D ane given by Eq. (3.123). This implies that, 
within the single-site KKR-CPA, the averaged Green’s function 
obtained by the augmented space formalism is exactly 
identical to that obtained by the familiar restricted 
averaging method. 

We now consider a two-atom cluster consisting of sites 
0 and 1. The configuration space, in this case, is spanned 

by four vectors, which are | f> = 1^1^ ~ ^^o ^1^’ 

if > = li.'° i->^> and !f...> = k? Tile augmented space for 

the cluster is spanned by eight vectors ji,/ i = 0, 1 and 
/ = f, f^, fj and f^j. For this case, Eq. (3.50) can be 
rewritten as 
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(c" - 


= 0,1 


i j=0 , 1 

/ ^ f 


i><i 1® \ 


(4.38) 

1=0,1 


f / V f 


, 1 J 

2 

\ oik 
Lu 

.|-(0,l)klQlj . . 

'eff ^eff 

(4.39) 


who re 


kl?=0, 1 

It IK evident that Also, it is evident 

from Eq. (4.38) that, ^ is a (6x6) matrix. The KKR-CCPA 

equations, given by Eq. (3.66) and (3.67), now reduce to the 

67 

following two equations: 


‘^eff = C - a> <0:fjrf;^|0;f^> 0) = c - u q““ u, 


= Q <0;f 03 = w 03, 

O ' iC ' i 


where 


00 -1 

Q = Rg 


(4.40) 

(4.41) 

(4.42) 


H- V- . K-> 


(4.43) 



Ill 


lh<* matricoK to are given by 


R = c - 
1 2 ’ 


R^ = C - R-^ 


„ I 4,00 ^-1 

Rj = C - *2 - “ K 2 


» r r>-l 01 

= C - 3^2 - V Rg V , 


R = C - ae°° - Oi r"^ Q, 
o 2 4 


(4.44) 


where 


H = ^ - I°° - R-l v'“, 
6 2 5 


.,^.3 ■P ‘ • n 1 

V = B *2 ’ ^ ’ 0 = 0 , 1 


(4.45) 


i j 

The variables Bfg can be found as follows. We partition 

the effective medium such that, the sites 0 and 1 fall inside 

the cluster. Now the effective medium path-operators in the 

67 

cluster subspace can be obtained by solving, 


' .^00 
'eff 

.j-01 ' 

‘ eff 


■ 

ef f 

- 1 
^eff 


[C 

*01 1 

*2 


.^10 
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Kq . (4.46) readily gives us, 


= ' 

2 

“eff ~ \} 

rOO _ 01 00 -1 

eff ~ eff efr 

. .flO 1-1 
ef fj ’ 

{4,47a) 

«01 

X = - 

2 

'^eff 

^00 01 -1 00 
'eff^'efr 'eff 

- 1"^ 
effj ’ 

(4.47b) 

«io 

*2 = - 

^eff * 

■-00 -10 .-1 00 
eff^'efr 'eff 

.^01 1-1 
' effJ ‘ 

(4.47c) 


We may point out that, the effective inverse t-matrices 
(C^f^) within the KKR-CCPA may not be diagonal any longer, if 
the cluster-size is beyond one-atom. Thus, care must be 
taken when evaluating the path-operators. 

4.4 RESULTS AND DISCUSSION 

In this section we present the results of our 
calculations of density of states for the one-dimensional 
muffin-tin alloy described in section 4.2. The lattice 
parameters (a = 6.00 a.u. ) and the muffin-tin radii (r^ = 
2.25 a.u. ) of the two components of the alloy are taken to be 
identical. The depth of the two constituent potentials are 
V, = -0.3 Ry and V„ = -0.5 Ry respectively. The densities of 
states for the two pure systems are shown in Fig. 4,1, which 
show peaks at band edges, characteristics of the 
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Fig. 4.1 Density of states for the one-dimensional 
chain of muffin-tin potentials of depth (a) 

0.3 Ry and (b) = 0.5 Ry. The lattice parameter 

B 

and muffin-tin radius are respectively 6.00 a.u 
and 2.25 a.u. for both the constituents. 
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<.ru«-dimcnKional model. For calculating the density of states 

within the KKH-CCPA, we have taken a two-atom cluster. We 

have calculated the densities of states for the A B alloy 

X y 

for t.wo ronc!entrat i ons (x = 0.1 and x = 0.5). 

We first discuss the low concentration limit. Fig. 
4.2(a) shows the density of states calculated within the 
KKR-CPA and KKR— CCPA. We note that, in the first majority 
band, there is no apparent difference in the KKR-CPA and 
KKR-CCPA resxilts. However, in the impurity band, there is a 
noticeable difference between the two. In contrast to a 
smooth KKR-CPA impurity band centered around E = -0.14 Ry, in 
the KKH-CCPA density of states we observe a shoulder in the 
KKR-CCPA density of states around E = -0.16 Ry and two sharp 
peaks at E = -0.12 Ry and E = -0.10 Ry respectively. The 
extra structures in the KKR-CCPA density of states, seem to 
ari.se due to correlated scattering from two-atom clusters. 
This becomes clear by a close inspection of the local density 
of state.s on impurities for a two-impurity cluster embedded 
in the effective medium. In Fig. 4.2(b), we show the local 
density of states on A-type impurities embedded in a pure 
B-medium. We observe that the two-impurity levels at E = 
-0.19 Ry and E = -0.10 Ry are close enough to the structures 
in the minority band of the KKR-CCPA density of states. The 
sharp structure at E = -0.12 Ry in the KKR-CCPA density of 
states is very close to the single-impurity level at E = 
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Fig. 4.2 (a) KKR~CPA (dashed line) and KKR-CCPA (solid line) 
density of states for the alloy with concentration x = 0.1. 
(b) Local density of states on the impurity site for a single 
(dashed line) and two (solid line) impurities of type A 
embedded in the pure B medium. 
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-0.13 Ky . W»- 111 so show tht* local density of states on 
imi»uritic's crabodded in the KKR-CPA medium in Figs. 4.3(a) and 
4.3(b). If we look at the local density of states on an 
A-type impurity of an AA-type clustert as shown in Fig. 
4.3(a), we observe that its impurity band is almost identical 
to that, of the KKR-CCPA density of states. The only 
difference is that in the KKR-CCPA density of states we have 
an extra structure at E « -0.12 Ry, which is very close to 
the peak at E = -0.13 Ry in the local density of states on 
the A-type impurity of an AB-type cluster, as shown in Fig. 
4.3(b). Also, a small but observable structure at E = 0.06 
Ry in the KKR-CCPA density of states is very close to the two 
kinks around E a 0.04 Ry, in the local density of states on 
both A-type and B-type impurities of an AB-type cluster, as 
shown in Fig. 4.3(b). The above discussion suggests that the 
structures in the KKR-CCPA density of states are due to the 
correlated scattering from two-atom clusters. 

We have also compared our KKR-CCPA results with some 
calculations obtained by using the embedded cluster method of 
Gonis et. We have embedded clusters of two, three, and 
five atoms in the KKR-CPA medium for x * 0.1. The averaged 
density of states on the central site of the cluster is 
obtained by calculating the local density of states on the 
central site for all possible configurations of the cluster 
and then taking the weighted average. This averaged density 
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(b) X =01 

LDOS on impurities 

On A 
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-0-4 0-0 0-4 0.^ 

Energy (Ry) 

Fig. 4.3 (h) Local density of states on the impurity site for 
a cluster of two impurities of AA (solid line) and BB (dashed 
line) type embedded in the KKR-CPA medium for x = 0.1. b) 

Local density of states on the impurity site of A < solid 
line) and B (dashed line) type of an AB-type cluster embedded 

in the KKR-CPA medium for x = 0.1. 
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of state:, [irtiiiui <•;; Konif* of the structures of the KKR-CCPA 
<ic>nsify ot states. In FiK.s. 4.4(a), 4.4(b), and 4.4(c), v?e 
show the averaged density of states, for different cluster 
.sizes. We observe t.hat., for a two-atom cluster, the averaged 
d<*nr. ity of states does not, produce the structures of the 
KKH-C-t'PA den.sity of states. However, as we increase the 
c 1 ust er-s 1 /.<* , the st ructures in the averaged density of 
states b<‘rom<* prominent, and show resemblance to the KKR-CCPA 
results. Ihi- .structure at E = -0.12 in the KKR-CCPA density 
of states is clearly reproduced when the cluster-size is 
incr<‘ase<l to five atoms. This further supports our assertion 
that the structures in the KKR-CCPA density of states are due 
to the c:orreInt.ed scattering from clusters of atoms. 

In Figs. 4.5(a) and 4.5(b), we show the component 
density of states and ^-decomposed KKR-CCPA density of states 
respectively. It is observed from Fig. 4.5(b) that, the 
lowe.st band is dominantly / = 0 in character, while the next 
higher band is / = 1 in character. 

In Fig. 4.6 we compare the KKR-CCPA density of states 
and the concentration-weighted average of the component 
density of states of Fig. 4.5(a). This weighted average is 
what one would get. by the restricted averaging method. It is 
important to note that these two densities of states are not 
same. It is expected, since the KKR-CCPA density of states 
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Fig. 4.4 Averaged density of states on the central site of a 
cluster of (a) two atoms, (b) three atoms, and (c) five 
atoms, embedded in the KKR-CPA medium for x = 0.1. 





Density of states { Arb. units) 



Energy (Ry) 

F i i; . 4.5 («) Component density of states on an A-type (solid 
line) and K-type (dashed line) atom in the alloy with 

concentration x = 0.1. (b) ^-decomposed density of states 

showing f = 0 (solid line) and ! =1 (dashed line) components 
of the averaged density of states for the alloy with 

concentration x - 0.1. 
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is obtained using the augmented space formalism by exactly 
averaging over the configurations of a two-atom cluster. 

Now we consider the case of a concentrated alloy (x = 
0.5). In Fig. 4.7, we show the densities of states 
calculated within the KKR-CPA and KKR-CCPA. In this case 
also, we observe that, the KKR-CPA gives a rather smooth 
density of states besides a peak in each band while the 
KKR-CCPA, once again gives few extra structures in all the 
bands. The KKR-CPA peaks at E = -0.39 Ry and E = -0.21 Ry 
almost coincide with the major KKR-CCPA peaks at E = -0.39 Ry 
and E = -0.21 Ry. Besides these major peaks, we observe some 
more peaks at E = -0.35 Ry, E = -0.29, E = -0.11 Ry, and E = 
-0.05 Ry, in the KKR-CCPA density of states. As in case of 
x=0.1, these extra structures in the KKR-CCPA density of 
states, are very close to the structures in the local density 
of states on the impurities of a two-impurity cluster 
embedded in the KKR-CPA medium. In Figs. 4.8(a) and 4.8(b), 
we show respectively the local density of states on the 
impurities of a cluster of two like and unlike impurities 
embedded in the KKR-CPA medium. We observe that, the 
KKR-CCPA peak at E = -0.35 Ry is exactly at the position of 
the peak in the local density of states on a B-type impurity 
of an AB-type cluster, as shown in Fig. 4.8(b). Thus the 
KKR-CCPA peak at E = -0.35 Ry may be identified as due to the 
B-type atoms. A small but observable structure at E = -0.29 
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Fig. 4.7 
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KKR-CPA (dashed line) and KKR-CCPA (solid line) 
of states for the alloy with concentration x = 0.5. 
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Ry in the KKR-CCPA density of states is at the exact location 
of a small peak in the local density of states on an A-type 
impurity of an AB-type cluster (Fig. 4.8(b)). The peak at E 
= -0.11 Ry is close to the peak in the local density of 
states on an A-type impurity of an AB-type cluster, which is 
at E = -0.14 Ry. One more extra peak at E = -0.05 Ry can be 
identified as due to the A-type atoms as it is exactly at the 
same location as a peak in the local density of states on an 
A-type impurity of an AA-type cluster, as shown in Fig. 
4.8(a) . 

In Fig. 4.9, we show the results corresponding to Fig. 
4.4, but for a different concentration, x = 0.5. As in the 

low concentration limit, in this case also, we observe that, 

the averaged density of states calculated by the embedded 
cluster method, does reproduce some of the structures of 
the KKR-CCPA density of states. The structures in the 
density of states of the embedded cluster method become 
prominent and sharper, and show an increasing resemblance to 
the KKR-CCPA density of states with the increase in the 

cluster-size. The extra structures in the KKR-CCPA density 

of states at E = -0.35, E = -0.29, and E = -0.11 are 
reproduced in the density of states of the embedded cluster 
method, when we take a five-atom cluster. Therefore, we 
conclude that, the structures in the KKR-CCPA density of 
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Fig. 4.9 Averaged density of states on the central site of a 
cluster of (a) two atoms, (b) three atoms, and (c) five 
atoms, embedded in the KKR-CPA medium for x = 0.5. 
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states are due to the correlated scattering from clusters of 
atoms . 

In Figs. 4.10(a) and 4.10(b), we show the component 
density of states and -f-decomposed KKR-CCPA density of states 
respectively for x = 0.5. In this case also, we observe that 
the lowest band is dominantly / = 0 in character. 

In Fig. 4.11 we compare the KKR-CCPA density of states 
and the concentration-weighted average of the component 
density of states presented in Fig. 4.10(a). Apart from the 
small structure at E = -0.29 Ry in the KKR-CCPA density of 
states, both these density of states seem to look alike. 
But, nevertheless, the height of the peaks are not same in 
the two cases. 

Just to see the effect of the different choices of 
off-diagonal terms between a real and an effective atom, we 
have calculated the local density of states on a single 
impurity embedded in the KKR-CCPA medium using Eq. (4.24), 
for various impurity potentials with different values of fJ . 
In Fig. 4.12(a), we show the local density of states on 
A-type impurity (potential-depth = 0.3 By) for x=0 . 1 , and p = 
0.0, 0.5, and 1.0. We observe that, there is some variation 
in the height of the peaks in the local density of states 
with fJ, although there is little qualitative change in the 
density of states. The variation in the height of the peaks 
for /J = 0.0 and fJ = 1.0 is less than 10%. A similar trend is 
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Fig. 4.10 (a) Component density of states on an A-type (solid 
line) and B-type (dashed line) atom in the alloy with 
concentration x = 0.5. (b) /-decomposed density of states 
showing / = 0 (solid line) and / =1 (dashed line) components 
of the averaged density of states for the alloy with 
concentration x = 0.5. 
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Fig. 4.11 KKR-CCPA density of states {solid line) and the 
concentration- weighted average of the component density of 
states (dashed line) for the alloy with concentration x = 

0 . 5 . 
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Fig. 4.12 Local density of states on a single impurity of 
potential -depth (a) 0.3 Ry (type A), and (b) 0.2 Ry, embedded 
in the KKR-CCPA medium for x = 0.1, using Eq. (4.24) for 
different values of u. 
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observed in the case of x = 0.5. Also, for an impurity of 

foreign kind, the same conclusion can be drawn from Fig. 

4.12(b), which shows local density of states on the impurity 

of potential-depth 0.2 Ry. It is, therefore, clear that the 

basic structure of the local density of states does not 

change with different choices of the off-diagonal terms 

67 

between real and effective atoms. 


4.5 CONCLUSIONS 


In this chapter, we presented the density of states for 
a one-dimensional muffin-tin alloy alloy calculated using the 
KKR-CCPA formulation presented in Chapter 3. We observed 
that the density of states calculated within the KKR-CCPA is 
more structured than the KKR-CPA density of states. We were 
able to correlate these structures with those in the local 
density of states on the impurities embedded in the KKR-CPA 
medium. We also found that these structures were missing in 
the averaged density of states calculated by the embedded 
cluster method for small clusters. However, when the 
cluster-size was increased, the embedded cluster method did 
reproduce some of the structures of the KKR-CCPA density of 
states. These observations suggest that, the extra 
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structures in the KKR-CCPA density of states arise due to 
correlated scattering from clusters of atoms in the effective 
medium. 

Our calculations of the local density of states on a 

single impurity embedded in the KKR-CCPA medium revealed that 

the basic structure of the local density of states remains 

unaltered with different choice of the off-diagonal terms 

6 7 

between real and effective atoms. However, we feel that, 
further work is needed to clarify the problem of embedding in 
the KKR-CCPA medium. 



CHAPTER 5 

DISORDERED ALLOYS WITH SHORT-RANGE ORDER 


5.1 INTRODUCTION 

The Korringa-Kohn-Rostoker cluster coherent-potential 

approximation (KKR-CCPA) presented in Chapter 3 is valid for 

6 7 

purely random alloys. Many disordered alloys are known to 

21 22 

exhibit short-range order. ’ By short-range order we mean 

the tendency of atoms of a given kind to surround atoms of 

the same kind or of a different kind. The augmented space 

formalism discussed in Chapter 2 , which was used to develop 

the KKR-CCPA formulation, has the restriction that, the site 

52 54 

occupation variables be statistically independent. ’ 

In other words, the probability of a particular kind of atom 

occupying a given site was assumed to be completely 

independent of the surrounding environment. However, for 

disordered systems with short-range order, this assumption is 

not valid. Disordered systems with short-range order can 

only be described in terms of dependent random variables. 

Thus, the augmented space formalism, presented in Chapter 2, 

can not be applied to these systems. However, the augmented 

space formalism can be generalized so that functions of 

, 70,94,95 

dependent random variables can also be averaged. 
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This new formalism has already been applied to disordered 

70 

alloys in the tight-binding framework. Our aim in this 

chapter is to extend this technique to the KKR-CCPA 

formulation presented in Chapter 3. 

We have applied the generalized augmented space 

formalism to obtain an expression for the configuration 

96 

averaged Green’s function for systems with short-range 
order. This averaged Green’s function is exact, and 

therefore, analytic. This is essentially an expansion of the 
configuration averaged Green’s function in presence of 

short-range order about that for a random system, and 

• 70,96 ^ , 

contains an infinite series. For numerical 

computations, we need to truncate this series after a few 

manageable terms. In the present work, we truncate the 

series after the first correction term. Using this 

formulation, we have calculated the density of states for the 

one-dimensional muffin-tin alloy with Markovian type 

70 

short-range order. We find that, our approximation yields 
non— negative density of states at all energies for a 
reasonably wide range of short-range order parameter. 
However, when the short-range order is strong, the density of 
states becomes negative in certain energy regions. This may 
be due to the approximation, in which we neglected the higher 
order correction terms in the expansion of the configuration 
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averaged Green’s function. The higher order terms may be 

important in the strong short-range order regime. We have 

compared our results with those obtained by the embedded 

50 

cluster method of Gonis and Freeman. We find that there is 

a general agreement between the two results. 

In section 5.2, we outline some mathematical 

preliminaries which we will need in the development of the 

theory. In section 5.3, we derive an expression for the 

configuration averaged Green’s function for disordered alloys 

with short-range order. In section 5.4, we describe the 
70 73 

Markov chain ’ and we present our results in section 5.5. 
Finally, in section 5.6, we outline the limitations and 
shortcomings of the theory. 


5.2 MATHEMATICAL PRELIMINARIES 

In Chapter 2, we mentioned that, a particular 
configuration of a disordered alloy can be described in terms 
of a set of occupation variables {n^}, in a particular 
sequence, and that a certain probability density is 
associated with each of these configurations. The simplest 
case that represents a disordered alloy is one in which {n^} 
are independent random variables. For this case, it is 
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possible t.o decompose Wn.D as a product of individual 

52 

probability densities of all the sites, i.e. 


= TT ( 5 . 1 ) 

1 

However, when (n^} are not statistically independent, 

such a decomposition is not possible. But, nevertheless, it 

is possible to decompose ip({n.}-) into two parts such that, 

^ 70 

the short-rani?e order part is separated. In the following 
we give the mathematical justification of the above 
statement . 

The probability density, dn. of a particular 

i 

configuration is a ’measure’ of how probable that particular 
configuration is, i.e. 


dfl = TT 

i 

97 

where {J is the measure. It is evident that, the measure 
will be different depending on whether {n^} are independent 
or dependent random variables. We define 


dpi = ID , ({n.OTT ’ 

^random i ' ! i 
1 

and 

du = iP ({n.D TT dn. , 

^2 1 


( 5 . 3 ) 


( 5 . 4 ) 
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where probability density in presence of 

short-range order. We note that, in a random distribution, 

all the configurations are ’probable’, however, in presence 

of short-range order, some or other configurations may be 

missing. In other words, if a configuration is allowed in 

the presence of short-range order, it must be present in a 

70 

random distribution. In a random distribution, the only 

configurations missing are those which are not consistent 

with the concentrations of the constituent atoms. Obviously, 

these configurations must be absent in presence of 

short-range order as well. In terms of measure, it can be 

stated as, if =0, then = 0. However, the reverse is 

not true. In measure theory, this condition is stated as, 

fJ^ is ’absolutely continuous’ with respect to . A basic 

theorem about the absolute continuous measures in measure 

97 

theory is the Radon-Nikodym theorem. When we apply this 

theorem to the present problem, it implies that it is 

70 

possible to write, 

<3/^2 = d^Jj^ (5.5) 

where ton.}) is called the Radon-Nikodym derivative of ^ 

^ .70 

with respect to jJ^ . Therefore, we may write, 
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i 

(5.6) 

The advantage of Eq. (5.6) is that, it allows one to 

calcuJate configuration averaged properties of a system with 

short-range order in terms of averages which are computed 

70 

with respect to independent random variables. The price 

one pays for this simplification, is the introduction of an 
unknown function (|3({n^)-) . However, we will show that, it 
does not present any difficulty. 

5,3 THE CONFIGURATION AVERAGED GREEN'S FUNCTION 


We consider a general function of dependent 

random variables (i.e. with short-range order). The 

configuration average of EG^^}) is given by, 

E«n D TT 

i 

For the configuration average of a function for a random 
system we will use the notation i which we have been using 
so far. 

Substituting Eq. (5.6) in Eq. (5.7), we get. 




Ea»i» TT Pilnjldn. ■ 


(5.8) 
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Now we are in a position to apply the augmented space 
52 

theorem given by Eq . (2.19) to Eq. (5.8), which gives us, 



= <f| Ea«S) |f>, 


(5.9) 


where ECiM^) and are obtained by replacing {n^} in 

EGn }) and (f<^n ^ by In order to use Eq. (5.9), we 

<X jL 

need to know I O- We note that, is a function 

of only, and therefore, is an operator in the 

configuration space alone (i.e. it is independent of the real 

Hilbert space). Thus, we can expand I f> in terms of 

the complete set of basis vectors in the configuration 
70 . 

space, i.e. 


where 




a„ = <f^l 


Substituting Eq. (5.10) in Eq. (5.9) we get, 


Hsro = 1 V ^ 1 


0a 


a 


(5.10) 


(5.11) 


(5.12) 
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The coefficients can be easily determined from known 
joint distribution functions. Putting ECiri^O = 1, in Eq. 

(5.12), we immediately get 

= 1 . (5.13) 

Putting ^5.12) we get 

<n.> = ^ <fl M" lf^>. (5.14) 

a 

Using Eq. (3.44a), we immediately get 
<n.> = X ag 4 a. , 

which readily gives us, 

a =0, for all i, (5.15) 

1 

because, Similarly, if we put = ^i^j 

Eq . (5.12), we get 

<7 


which with the help of Eqs. (3*44a) and (3#44c)i gives us, 
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2 1/2 

<n.n.> = a. x + x (xy) (a. + a.) + xy a. . . 

1 .1 ^ 1 J 1 J 


(5.16) 


With the help of Eqs. (5.13) and (5.15) we get, 


a., = [<n.n.>-x ] /xy , 

1 j 1 j 


(5.17) 


From the definition of {n.} (Eq. (2.32)), we get 

1 


<n . n .> = , 

a .1 AA 


(5.18) 


i.e. it is the probability of having A-type atoms on i-th and 


i j 2 

j-th sites. In a random distribution = x , so that a^ 


AA 

vanishes in this case. In the present case, we have 




a. . = [ P^^ - ] /xy. 

1 J AA 


(5.19) 


In a similar manner, we can evaluate all the 
7 0 

coefficients. In the following we list a few of them. 


r 

a... = LP... " X 

1 J k AAA 


AA AA AA 

+ CP 

2 


i jkl *■ ' AAAA 

^jk _kl ^li 
AA 

and so on 


(5.20) 


. _ijkl X X + X^ fp^’^ 

a. = [Py,,jAA - X CP^^^ AAA AAA AA 


P - . p*-^ . p“ H . P^^ . pfb -3/ ]/(xrt 

AA AA AA AA AA 


(5.21) 
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The: configuration averaged Green’s function can be 

written from Eq . (5.12) as, 

[<G>]gj^O (5.22) 

a 

This is an infinite sum. The first term in the series, 

which is ^fjGv'i^ }')|0» is the configuration averaged Green’s 

function for a purely random system. Therefore, Eq. (5.22) 

can be regarded as a perturbation expansion, when the higher 

terms in the series come as a perturbation to the averaged 

Green’s function of the random distribution. Thus, Eq. 

(5.22) allows us to calculate the averaged Green’s function 

in presence of short-range order, in terms of the random 
70 

averages . 

If we apply this technique to the site-diagonal form of 
the KKR Green’s function, we will very easily get. 


SRO 


<6sp> 


+ Tr 


C^0 


B 


<o:f| 'Sf]. <5-23) 

where <6„„> is given by Eq. (3.116). Since Eq. (5.23) 
involves an infinite sum, it can not be used for numerical 
computations. Therefore, we have to introduce some 
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approximation, i.e. to terminate the sum after a few steps. 
In our calculation, we have terminated the sum after two 
steps. In this case, only a. . are non-zero (because, all 

X 0 

a^’s are zero, Eq . (5.15)). When we apply the CCPA to the 
present case, the number of terms to be summed can be further 
reduced. For a two-atom cluster consisting of sites 0 and 1, 
the only non-zt'ro term in the expansion corresponds to . 
This essentially means that, we keep only the first 
correction tei'm in the expansion of the configuration 
averaged Green’s function about the averaged Green’s function 
of the random system. With this approximation, Eq. (5.23) 
can be written as 


<Gs 0 > 


= <^G y + 

SRO 


K - •■] " [«' 


.-1 


-B 


flA 10;fQ^>F‘' + 


,-l 

<0;flA M°l0;fQ^>(5F 




(5.24) 


Using Eq. (3.44d), we get 


< Gsd > 


SRO 


<G > + -^ - x^l Tr [<0;f|A 

xy [ AA J L 


^-1 


<0 


;f|A |o;fj>6Fj, 


0;f„i> F * 


(5.25) 


Z. ,-A i-B 
F = y F + X F • 


where 
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The matrix elements of A ^ appearing in Eq. (5.25) can 
be calculated easily, 


✓a.^'Ia Ia.^- X ft' 00 , 01 , OU 01 

^0,f]A I ’^01^ ~ ^^eff ~ ^2 ^ ^ “ (V + b ) oo Q , 


(5.26) 


and 


-1 


<0;f|A 10;fj> = 


-01 , 01 . .10 
(v + b ) 0 ) b 


(5.27) 


The variables appearing in Eqs . (5.26) and (5.27) are 

defined in Chapter 4 (section 4.3). Therefore, the only 

unknown quantity in Eq. (5.25) is , which we will define 

AA 

in the next section. 


5.4 MARKOV CHAIN 

This is the simplest possible system, for which 

short-range order can be easily defined. The name comes from 

a well-known process in the probability theory, called Markov 

process In a Markov process, the probability of the i-th 

event is restricted to depend only on the probability of the 

70 

(i-l)-th event. By a close analogy, in a Markov chain, the 
occupation probability of i-th site is assumed to depend only 
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on t.ho kind cjf atom occ:upying its first nearest— neighbor 

sites. Therefore, in a Markov chain, the short-range order 

is restricted to nearest -nei ghbor pairs only. Following 
. , 4 6 

Kumar e£ a! , we define the nearest-neighbor pair 
distribution as: 

P<n. , n. = 1) = (1- ay) 6(n. - 1) +ay6(n.), (5.28) 

d ^ 0 J 

and 

p(n. , n. = 0) =ax<5(n. - 1) + (1- ax) 6(n.). (5.29) 

^ 0 3 

Here a is the parameter that describes the short-range 
order. It is easily seen that, a = 1 represents a system 
with no short-range order, for a < 1, AA- and BB-type 
clustering is favored, and for a > 1, AB-type ordering would 
be dominant. The maximum value a = 1/y corresponds to a 

situation, when the alloy is not random at all, but forms a 
superlattice of A and B constituents. The minimum value a = 
0 corresponds to a situation, when A- and B— type atoms are 
totally segregated. When a is close to unity, 
short-range order is weak. 

From Eq. (5.28), we easily get 

01 , ^ 

P =x(l-ay)=xy(l-a)+x. 

AA 


(5.30) 
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We wiJJ use this form of P in Eq. (5.25), for the 
calculation of the density of states of this system. 


5.5 RESULTS AND DISCUSSION 

We have calculated the density of states for the one- 
dimensional muffin-tin alloy with Markovian type short-range 
order. In Fig. 5.1 and 5,2, we present the results for x = 
0.1 for various values of short-range order parameter oi. 
Fig. 5.1(a) shows the density of states for Cf. = 1.0, which 
corresponds to the random alloy. In Fig. 5.1(b), we show the 
density of states for a = 1.11, which is the maximum allowed 
value of Ct for X = 0 , 1 . We observe that, there is little 
difference between this density of states and the density of 
states for a random system. In Figs. 5.2(a) and 5.2(b), we 
show the density of states for a = 0.5, ot = 0.75, and a = 
0.25 respectively. We note that, for a = 0.25, the density 
of states becomes negative in the region between E = -0.17 Ry 
and E = -0.09 Ry. The range of a, for which the density of 
states is non-negative throughout the whole energy range is 
from 0.25 to 1.11, which is appreciable. We observe a 
general trend in these results. We note that, when a 
decreases gradually, some states are pushed into the first 
majority band, while the impurity band goes on decreasing 



Density of states (Arb. units) 



Energy (Ry) 

Fig. 5.1 Averaged density of states for the alloy with short- 
range order, for the concentration x = 0.1. The short-range 
order parameters are (a) oi = 1.0 and (b) ot = 1.11. 
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(a) X = 01 

c^=0-5 

(/ = 0.75 





(b) X =01 .c/ = 0-25 


“04 


00 

Energy (Ry) 


Fig. 5.2 Averaged density of states for the alloy with short- 
range order, for the concentration x = 0.1. The short-range 
order parameters are (a) « = 0.75 (dashed line), ot - 0.50 
(solid line), and (b) 0( = 0.25. 
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with 0( . This is because, when 0( decreases, the probabilities 

of BB and AA pairs increases, while that of an AB pair 

decreases. For x = 0.1, the probability of a BB pair is 

quite large compared to the other clusters. Hence, with the 

decrease in a, the system tends towards a pure B system. 

Therefore, the impurity band, which is the characteristics of 

A-type atoms, decreases gradually with the decrease in O'. We 

have compared our results with those calculated by the 

50 . 

embedded cluster method of Gonis and Freeman, in which one 
can include the effects of short-range order, but in a 
non-self-consistent way. In Figs. 5.3 and 5.4, we show the 
density of states obtained by this method, for various values 
of a. In this case also, we observe that, with the decrease 
in a, the impurity band loses some states, while the first 
majority band gains. Thus, there is a general agreement 
between the two results. 

In Figs. 5.5 and 5.6, we show the density of states 
calculated by our formulation for x = 0.5. In Figs. 5.5(a), 
5.5(b), and 5.5(c), we show the density of states for O. = 1.0 
(random), a = 1.25, and O'- = 1.5 respectively. In this case 
also, we observe that, when a = 1.5, which is much away from 
the value corresponding to the random distribution (ot = 1.0), 
the density of states becomes negative in some energy 
regions. In Figs. 5.6(a) and 5.6(b), we show the results for 
cx = 0.5, a = 0.75, and a = 0.25 respectively. Again, we 


Density of states (Arb. units) 



(a) X =01 ,o^ = 1-0 



00 

Energy (Ry) 


Fig. 5.3 Averaged density of states for the alloy with short- 
range order, for the concentration x = 0.1 by the embedded 
cluster method of Ref. 50. The short-range order parameters 
are* (a) oi = 1,0 and (b) a = 1.11. 




Densit 



Energy (Ry ) 


Fig. 5.4 Averaged density of states for the alloy with short- 
range* order, for the concentration x = 0.1 by the embedded 
clu.ster method of Ref. 50. The short-range order parameters 
are (a) o - 0.75 (dashed line), a = 0.50 (solid line), and 






Density of states (Arb. unites) 


(a) X = 0 5 , o'M O 



(b) X = 0 5 , o<'=l-25 



Fig. 5.5 Averaged density of states for the alloy with sho 
range order, for the concentration x = 0.5. The short-ra 
order parameters are (a) ct = 1.0. (b) a = 1.25, and (c) ct 
1.50. 
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observe that, for 0i = 0.25, the density of states becomes 
noKative in some energy regions. In this case also, the 
range of a, for which we get non-negative density of states 
is appreciable (from 0.25 to 1.50). We find that, as cx 
decreases, states are pushed into the first band. This is 
expected, since the first band is the characteristics of 
B-type atoms, and with the decrease in O', the probability of 
BB pairs increases. In Figs. 5-7 and 5.8, we show the 
corresponding results calculated by the embedded cluster 
method. We note that, there is a general agreement between 
these results and the results obtained by our formulation. 

It is clear from the above discussion that the 
approximation in which we have neglected the higher order 
correction terms in the expansion of the configuration 
averaged Green’s function, gives non-negative density of 
states for a considerably wide range of the short-range order 
parameter a. It gives negative density of states only when oi 
is very far from unity, when the higher order terms in the 
expansion become important. It may be mentioned that, one- 
dimensional systems are known to be very sensitive to 
approximations. The fact that, our approximation is 
reasonably good for a on£‘-damensional system, gives one hope 
that., it will be quite satisfactory when applied to a three- 
dimcTjsional system. 



X = 0-5 , 0(^ = 1 00 


(b) X = 0-5 , 0^ = 125 


(c) X =0-5 , o^ = 1-50 


-0 4 00 0 4 01 

Energy (Ry) 

Fig. 5.7 Averaged density of states for the alloy with short- 
range order, for the concentration x = 0.5 by the embedded 
f m(*thod of Ref . 50. The short-range order parameters 

are (a) = 1.0, (b) c* = 1.25, and (c) a = 1.50. 










5.6 CONCLUSION 


In this chapter, we presented the generalization of the 
configuration averaging process in the augmented space 
formalism for the systems with short-range order. We used 
this formalism to determine the configuration average of the 
site-diagonal form of the KKR Green’s function. This 
averaged Green’s function is exact, and therefore, is 
herglotz. However, for numerical computations, one has to 
adopt some approximation. This approximation is in addition 
to the CCPA, which is used to reduce the dimensionality of 
th<‘ augmented space. We calculated the density of states for 
a one-dimensional muffin-tin alloy with Markovian type 
short-range order. We found that this approximation yields 
non-negative density of states for a considerably wide range 
of values of the short-range order parameter. 



CHAPTER 6 

SUMMARY AND CONCLUSIONS 


In this thesis, we discussed a self-consistent cluster 

3-5 

generalization of the coherent-potential approximation 

(CPA) for disordered substitutional alloys. This cluster 

generalization, which is called the cluster coherent- 

46,58 

potential approximation ( CCPA ) , is developed in the 

52-55 

augmented space formalism. This approximation includes 

the correlated scattering from atoms inside a chosen cluster, 
in a self-consistent way, which may be important in systems 
showing short-range order and clustering tendencies, such as 
CuPd , CuPt, and CuNi systems. The augmented space CCPA 
restores the translational symmetry of the effective medium. 
In this thesis, we proved three theorems regarding the 
analyticity of the augmented space CCPA. It was shown that, 
t hc‘ augmented space CCPA preserves the analytic (herglotz) 

pioperties of the self-energy and the configuration averaged 

5 S 

Green’s function. Also, it was proved that the 

self-consistent iteration scheme for the self-energy has a 

fixed point,, i.e. it always converges to a herglotz, unique 

and liounded value, provided we start off the iteration with a 

5 8 

tK)unded and herglotz one. This guarantees a single-valued, 
continuous, and non-negative density of states at all 


energ i es . 
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'I ho main objcctave of this thesis was to present a 

cluster generalization of the first-principles parameter-free 

67 88 

Korringa-Kohn-Rostoker coherent-potential approximation ’ 

{ KKR-CPA ) . We have developed a cluster theory, the KKR-CCPA, 

by effectively combining the augmented space CCPA and the 

conventional KKR Green’s function formulation of the 

electronic structure of disordered alloys. The KKR-CCPA 

formulation has been put on the same footing as the 

single-site KKR-CPA, in which one can calculate the charge 

densities, and therefore, can achieve charge self-consistency 

within the local density approximation of the density 

13-15 

functional theory. Therefore, our KKR-CCPA formulation 

can be regarded as a first-principles parameter-free theory 
of disordered substitutional alloys, in which correlated 
scattering from clusters of atoms is taken into account. The 
KKR-CCPA restores the translational symmetry of the effective 
medium, and like the KKR-CPA, is exact in the limit of 
vanishing concentration of either constituent. Also, it 

rt'duces to the correct limit (KKR-CPA), when we take a 

, 67 

one-atom cluster. 

The KKR-CCPA formulation was applied to calculate the 

density of states for a one-dimensional muffin-tin alloy. We 

fo\ind that, the density of states calculated within the 

6 7 

KKR-CCPA is more structured compared to the density of 
states calculated within the single-site KKR-CPA. To 

understand and investigate the origin of these structures, we 
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calculated the local density of states on an impurity of 
different two-impurity clusters embedded in the KKR-CPA 
medium. We found that, the structures in the KKR-CCPA 
density of states and those in the local density of states on 
these impurities were correlated. Also, we compared our 
results with the averaged density of states of the embedded 
cluster method of Gonis and coworkers, for various 

cluster-sizes. It was found that, with the increase in the 

cluster-size the averaged density of states obtained by the 

embedded cluster method showed more resemblance to the 

KKR-CCPA density of states. From these observations, it was 
concluded that, the structures in the KKR-CCPA density of 

states originate from the correlated scattering from clusters 

r . 67 

of atoms. 

We have formally extended the KKR-CCPA formulation to 

include non-random effects like the short-range order. We 

showed that the configuration averaged Green’s function in 

presence of short-range order, can be expanded about that of 

a purely random system. The correction terms could be 

expressed in terms of the averaged inverse t-matrices and the 

configuration averaged Green’s function for a purely random 
70 

system. For numerical computations, we adopted an 

approximation, in which we truncated the infinite series in 
t iie expansion after the first correction term. We applied 
this formulation to calculate the density of states for the 
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one-dimensional muffin-tin alloy with Markovian type 

short-range order. We found that, the density of states 

remained non-negative throughout the whole energy range, for 

an appreciably wide range of the values of the short-range 

order parameter. The density of states became negative at 

some points, only in the strong short-range order regime. 

This unphysical behavior is probably due to the breakdown of 

the approximation, in which the higher order corrections in 

the expansion of the averaged Green’s function are neglected. 

The contribution of the higher order terms may be significant 

in the strong short-range order regime. In this regime, the 

system is closer to an ordered system, rather than a 

disordered alloy. So, when the alloy is close to an ordered 

system, it may be more appropriate to expand the averagad 

99 

Green’s function about its value for an ordered system, 
rat.her than that of a random system. However, we found that, 
our approach is successful over an appreciably wide range of 
short-range order parameters for a one-dimensional model, 
which is known to be sensitive to any kind of approximation. 
Therefore, we hope that, our approach will be more successful 
for a three-dimensional system. 

In this work, we have applied the KKR-CCPA formulation 
to a model case. In principle, this theory can be applied to 
realistic alloy systems in three dimensions. However, 
comijutat ional 1 y it will be much involved than the KKR-CPA. 



162 


fkjt , rif ve rthe loss , we fee] that, with the availability of the 
modern computers, the application of the KKR-CCPA to 
realistic alloy systems can be realized. This problem has 
been taken up by our group, and we are quite hopeful that, we 
will be able to report the results within a reasonable time. 

We observed sharp structures in the density of states of 
the one-dimensional muffin-tin alloy calculated within the 
KKR-CCPA. We expect that, in three dimensions also, the 
KKR-CCPA will produce sharp structures in the density of 
states. The appearance of the structures in the averaged 
density of states may strongly influence some of the physical 
properties of the alloy. For example, the Fermi energy and 
the density of states at the Fermi energy are expected to be 
influenced by the appearance of sharp structures in the 
density of states, especially, if the structures are in the 
Vicinity of the Fermi energy. Any change in the Fermi energy 
or the density of states at the Fermi energy, will affect the 
magnetic and other properties of the alloy. The formation of 
magnetic moments, the electronic specific heat, 
superconducting transition temperature, electrical 
conductivity, Fermi surface etc. are among such properties. 
T]je above considerations underline the importance of the 
KKR-CCPA formulation. 

A very important and long sought goal in the theory of 

100 


disordered alloys is the calculation of phase diagrams 
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from first-principles. In the determination of phase 

diagrams by statistical methods, effective pair (and, if 

necessary, larger cluster) interactions are very 
101 

important. The effective pair interactions have been so 

far calculated by embedding two-atom clusters in the KKR-CPA 
102 

medium, which is a non-self-consistent approach. A better 

approach would be to use self-consistently calculated 

effective pair interactions, which come out naturally from 

our KKR-CCPA approach. Therefore, it is hoped that, our 

formulation will be helpful in accurate determination of 

phase diagrams from the first-principles. 

The application of the KKR-CCPA formulation to the 

systems with the tendency to clustering, will be more 

int.eresting , especially, in the investigation of those 

properties which depend intrinsically upon the clustering 

effects. In CuNi alloy, it is well-known that, Ni atoms will 

carry a magnetic moment only if surrounded by a sufficiently 

16 17 

large number of other Ni atoms. ’ In CuPd system, the 

dtujsity of states calculated within the KKR-CPA (Refs. 23 and 

26 

24) does not agree with the experimental results, 

27 

especially in the lower energy region. Many systems like 
CuPt and CuPd are known to have short-range ordering. These 
systems are, therefore, ideal systems for the application of 


the KKR-CCPA. 
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One* problem which needs further work, is the embedding 

of a cluster of real atoms in the KKR-CCPA effective medium. 

The KKR-CCPA, quite predictably, introduces site-off-diagonal 

corrections as well in the inverse effective t-matrices, 

along with the usual site-diagonal corrections. The 

off-diagonal corrections in the t-matrices are quite 

complicated and their dependence on the potentials of the two 

involved sites is not explicitly known. This poses problem 

in determining the off-diagonal term between a real atom and 

an effective atom, when a cluster of real impurities is 

embedded in the effective medium. Though, we found that, the 

basic structure of the local density of states on the 

impurities is almost independent of different choices of the 

6 7 

off-diagonal terms, we feel that, this problem requires 
some attention. 


The study of topologically (structurally) disordered 

alloys is largely untouched. The KKR-CCPA formulation 

presented in the present thesis also, is valid only for 

substitutionally disordered alloys. In most of the systems, 

especially in liquid alloys and amorphous solids, topological 

103 

disorder is coupled with strong short-range order. Some 
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lead to a herKiotz Green’s function, and analyticity 

difficulties and convergence problems show up in 
1 10 

calculations. We have tackled the problem of short-range 

order to some extent in this thesis. It will be interesting 
to generalize this formulation for systems with topological 
disorder. 


Our KKR-CCPA formulation is a non-relativistic theory. 

For disordered alloys containing heavy elements (like 5d 

elements), the relativistic effects are expected to be 

important. A fully relativistic version of the single-site 

KKR-CPA and its successful application to many systems have 

87 


been reported in the literature 


Very recently, it even 


became possible to develop a cluster generalization of the 

111 


relativistic KKR-CPA in the embedded cluster method. 


The 


relativistic KKR-CPA embedded cluster method, in essence, is 

. , 87 

a combination of the relativistic KKR-CPA of Staunton et at 

50,60 

and the embedded cluster method of Gonis and coworkers. 

This has been very successful in predicting the short-range 
order effects in CuPt (Ref. 112) systems. But, nevertheless, 
this remains a non-self-consistent approach. A 


self-consistent cluster generalization of the relativistic 
KKR-CPA may be realized, if one can efficiently apply the 
augmented space formalism to this problem. 


REFERENCES 


1. L. Nordheim, Ann. Phys. (Leipzig) 9, 607 and 641 (1931); 

T. Muto, Sci . Pap. Inst. Phys. Chem. Res. (Jpn.), 34, 377 

(1938); F. Bassani and D. Brust, Phys. Rev. 131 , 1524 

(1963); H. Amar , K. H. Johnston, and C. B. Sommers, Phys. 
Rev. 153, 655 (1967); M. M. Pant and S. K. Joshi, Phys. 
Rev. 184, 635 (1969). 

2. J. Korringa, J. Phys. Chem. Solids, 7, 252 (1958). 

3. P. Seven, Phys. Rev. 156, 809 (1967). 

4. D. W. Taylor, Phys. Rev. 156, 1017 (1967). 

5. B. Velicky, S. Kirkpatrick, and H. Ehrenreich, Phys. Rev. 
175, 747 (1968). 

6. H. Ehrenreich and L. M. Schwartz, in Solid State Physics, 
edited by H. Ehrenreich, F. Seitz, and D. Turnbull 
(Academic, New York, 1976), Vol . 31, p.l49. 

7. J. S. Faulkner, in Progress in Materials Science , edited 
by J. W. Christian, P. Hassen, and T. Massalski (Pergamon, 
New York, 1982), Vol. 27, p. 1. 

8. A. Bansil, in Electronic Band Structure and its 

Applications , edited by M. Yussouff ( Springer-Verlag , 
Berlin, 1987), p. 273. 

9. E, Muller-Hartmann, Solid State Commun. 12, 1269 (1973). 

10. F. Ducastelle, J. Phys. C 7, 1795 (1974). 

11. B. Velicky, Phys. Rev. 184, 614 (1969). 

12. J. S. Faulkner and G. M. Stocks, Phys. Rev. B 21, 3222 

( 1980) . 

13. W. Kohn and L. J. Sham, Phys. Rev. 140, A 1113 (1965). 

14. F* . Hohenberg and W. Kohn, Phys. Rev. 136, B 864 (1964 ). 

15. S. Lundqvist and N. H. March, in Theory of the 

Inhomogeneous Electron Gas, edited by S. Lundqvist and N. 
H. March (Plenum, New York, 1983); J. Callaway and N. H. 
March, in Solid State Physics, edited by H. Ehrenreich, F. 
Seitz, and D. Turnbull (Academic, New York, 1984), Vol. 
38 , p . 135. 



167 


H) . J. F. Perrier, B. Tissier, and R. Tournier, Phys . Rev, 
Lett. 24, 313 (1970). 

17. J, S. Kouvel and J. B. Comly, Phys. Rev. Lett. 24, 598 

{ 1970) . 


18. J. Inoue and M. Shimizu, J. Phys. Soc. 

Jpn . 

4, 1321 

( 1976) . 




19. J. C. Ododo and B. R. Coles, J. Phys. F 7, 

2393 

( 1977). 


20. K. J. Kim and D. W. Lynch, Phys. Rev. B 39 

, 9882 

( 1989 ) 

• 

21. J. Kulik, D. Gratias, and D. de Fontaine, 

Phys . 

Rev . 

B 

40, 8607 (1989). 




22. J. Banhart , W. Pfeiler and J. Voitlander, 

Phys . 

Rev . 

B 


37, 6027 (1988). 

23. R. S. Rao, A. Bansil, H. Asonen, and M. Pessa, Phys. Rev. 


B 29, 1713 (1984 ) . 

24. H. Winter, P. J. Durham, W. M. Temmerman, and G. M. 
Stocks, Phys. Rev. B 33, 2370 (1986). 

25. D. van der Marel, J. A. Jullianus and G. A. Sawatzky, 
Phys. Rev. B 32, 6331 (1985). 

26. H. Wright, P, Weightman, P. T. Andrews, W. Folkerts, C. 
F. J. FI ipse, G. A. Sawatzky, D, Norman and H. Padmore , 
Phys. Rev. B 35 , 519 (1987). 

27. N. Stefanou, R. Zeller, and P. H. Dederichs, Solid State 
Commun. 62, 735 (1987). 

28. M. Tsukada, J. Phys. Soc . Jpn. 27, 684 (1969); ibid., 32, 
1475 (1972). 

29. B. G. Nickel and J. A. Krumhansl, Phys. Rev. B 4, 4354 

( 1971 ) . 

30. F. Cyrot-Lackman and F. Ducastelle, Phys. Rev. Lett. 27, 
429 (1971 ) . 

31. K. F. Freed and M. H. Cohen, Phys. Rev. B 3, 3400 (1971). 

32. R. Mills, Phys. Rev. B 8, 3650 (1973). 

33. F. Brouers, F. Ducastelle, and J. van der Rest, J. Phys. 
F 3 , 1704 ( 1973 ) . 

34. F. Brouers and F. Ducastelle, J. Phys. F5, 45 (1975). 

35. N. Zaman and R. J. Jacobs, J. Phys. F 5, 1677 (1975). 


168 


36. N. F. Berk, D. J. Shazeer, and R. A. Tahir-Kheli, Phys. 
Rev. B 8, 2496 ( 1973) . 

37. N. F. Berk and R. A. Tahir-Kheli, Physica 67, 501 (1973). 

38. W. H. Butler, Phys. Lett. A 39, 203 (1972). 

39. V. Capek, Phys. Status Solidi B 52, 399 (1972). 

40. F. Brouers, F. Ducastelle, F. Gautier, and J. van der 
Rest, J. Phys, F3, 2120 (1973). 

41. F. Brouers, M. Cyrot, and F. Cyrot-Lackman , Phys. Rev. B 
7, 4370 (1973). 

42. K. Aoi, Solid State Commun, 14, 929 (1974). 

43. W. H. Butler, Phys, Rev. B 8, 4499 (1973). 

44. A. R. Bishop and A. Mookerjee, J. Phys. C 7, 2165 (1973). 

45. R. Mills and P. Ratanavararaksa , Phys. Rev. B 18, 5291 

( 1978 ) . 

46. V. Kumar, A. Mookerjee, and V. K. Srivastava, J. Phys. C 
15, 1939 (1982); P. K. Thakur, A. Mookerjee, and V. A. 
Singh, J, Phys. Fl7, 1523 (1987). 

47. A. Gonis and J. W. Garland, Phys. Rev, B 16, 2424 (1977). 

48. C. W. Myles and J. D. Dow, Phys. Rev. Lett. 42, 254 

( 1979) . 

49. C. W. Myles and J. D. Dow, Phys. Rev. B 19, 4939 (1979). 

50. A. Gonis and A. J. Freeman, Phys. Rev. B 29, 4277 (1984). 

51. B. G. Nickel and W. H. Butler, Phys. Rev. Lett. 30, 373 

( 1973 ) . 

52. A. Mookerjee, J. Phys. C6, L205 (1973); 6, 1340 (1973). 

53. T. Kaplan and L. J, Gray, Phys. Rev. B 14, 3462 (1976); 

ibid., B 15, 6005 (1977); J. Phys. C 9, L303 (1976). 

54. A. Mookerjee, in Electronic Band Structure and its 
Applications , edited by M. Yussouff ( Springer-Verlag , 
Berlin, 1987), p. 248. 

55. P. L. Leath, in Excitations in Disordered Systems, edited 
by M . F. Thorpe (Plenum, New York, 1982) p. 109; T. Kaplan 
and L. J. Gray, ibid., p. 129. 

56. A. Mookerjee, J. Phys. C9, 1225 (1976); 19, 193 
19, 275 (1986). 


( 1986) ; 



169 


57. T. Kaplan, P. L. Leath, L. J. Gray, and H. W. Diehl, 
Phys. Rev. B 21 , 4230 (1980). 

58. S. S. A. Razee, A. Mookerjee, and R. Prasad, J. Phys. 
Condensed Matter (in Press). 

59. A. Gonis, W- H. Butler, and G. M. Stocks, Phys. Rev. 
Lett. 50, 1482 (1983). 

60. A. Gonis, G. M. Stocks, W. H. Butler, and H. Winter, 
Phys. Rev. B 29 , 555 (1984). 

61. R. Mills, L. J. Gray, and T. Kaplan, Phys. Rev. B 27, 
3252 (1983). 

62. T. Kaplan and L. J. Gray, Phys. Rev. B 29, 3684 (1984). 

63. A. Mookerjee, Disordered Systems (Hindustan Publishing, 
Delhi, 1979). 

64. A. Mookerjee, J. Phys. F 17, 1511 (1987). 

65. S. S. Rajput, S. S. A. Razee, R. Prasad, and A. 

Mookerjee, J. Phys. Condensed Matter 2, 2653 (1990). 

66. A. Gonis and G. M. Stocks, Phys. Rev. B 25, 659 (1982). 

67. S. S. A. Razee, S. S. Rajput, R. Prasad, and A. 

Mookerjee, Phys. Rev. B 42, 9391 (1990). 

68. A. Mookerjee and P. K. Thakur , Phys. Rev. B 38, 3798 

(1988) . 

69. W. H. Butler, Phys. Rev. Bl4, 468 (1976). 

70. L. J. Gray and T. Kaplan, Phys. Rev. B 24, 1872 (1981). 

71. R. Haydock, in Excitations in Disordered Systems, edited 
by M . F. Thorpe (Plenum, New York, 1982), p. 29. 

72. J. A. Shohat and J. D. Tamarkin, The Problem of Moments 
(American Mathematical Society, 1963). 

73. D. Paquet and P. Leroux-Hugon, Phys. Rev. B 29, 593 

( 1984 ) . 

74. A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski , 
Methods of Quantum Field Theory in Statistical Physics 
(Prentice-Hall, New Jersey, 1963). 

75. P. R. Halmos, Hilbert Space Problem. Book (Van Nostrand, 
New York , 1967 ) . 

76. A. Bansil, L. M. Schwartz, and H. Ehrenreich, Phys. Rev. 
B 12, 2893 (1975). 



170 


77. R. Pfasad , S. C. Paf)adopou 1 os , and A. Bansil, Phys. Rev. 
B 23 , 2607 (1981). 

78. R. Pra.sad and A. Bansil, Phys. Rev. Lett. 48, 113 (1982). 

79. H. .^sonen , M. Lindroos, M. Pessa, R. Prasad, R. S. Rao , 
and A. Bansil, F’hys. Rev. B 25, 7075 (1982). 

80. R. Prasad, R. S. Rao, and A- Bansil, in Excitations in 

Disordered sys terras , edited by M. F. Thorpe (Plenum, New 

York, 1982), p. 241. 

81. R. Benedek, R. Prasad, S. Manrnnen, B. K. Sharma, A. 

Bansil, and P. E. Mijnarends, Phys. Rev. B 32, 7650 

(1985). 

82. R. S. Rao, R. Prasad, and A. Bansil, Phys. Rev. B 28, 
5762 (1983). 

83. R. Prasad, in Electronic Band Sirv.ctv.re and its 

Applications , edited by M. Yussouff ( Springer-Verlag , 

Berlin, 1987), p. 299. 

84. P. Lloyd, Proc. Phys. Soc . , London 90, 207 (1967). 

85. R. Prasad and A. Bansil. Phys. Rev. B 21, 496 (1980). 

86. S. Kaprzyk and P. E. Mijnarends, J. Phys. C 19, 1283 
(1986 ) . 

87. J. B. Staunton, B. L. GyOrffy, and P. Weinberger, J. 
Phys. F 10, 2665 (1980); P. Weinberger, A. M. Boring, R. 
C. Albers, and W. M. Temmerman, Phys. Rev. B 38, 5357 
(1988); E. S. Aerola, C. J. Barnes, R. S. Rao, and A. 
Bansil, Phys. Rev. B (to be published). 

88. S. S. A. Razee and R. Prasad (unpublished). 

89. In Ref. 67, the joint averages were evaluated by the 
restricted averaging method of Ref. 12. The restricted 
averaging method is not appropriate for the KKR-CCPA, 
because, in this method, the statistical correlation 
between the sites is not included, and therefore, it is 
essentially a single-site appro.ximat ion (Ref. 12). 

90. P. J. Durham, B. L. GyOrffy, and A. J. Pindor, J. Phys. F 
10, 661 (1980). 

91. W. M. Temmerman, J. Phys. F 12, L25 (1982). 




171 


92. 'I’C) determine t,he component charge densities, we embed a 

single atom in the effective medium. In principle, one 

can embed a cluster of impurities to find the charge 

densities for different configurations of the cluster. 

a -+ 

The component charge density p (r ) on an c-type atom 

will be different for different configurations of atoms 

around it. This, in turn, will give us different values 

of V (r ) for different configurations of the cluster. 

o 

•4 

However, in our formulation, we assume that o (r ) is same 

o 

for different clusters. 

93. E. Merzbacher, Quanturi Mechanics (Wiley, New York, 1970). 

94. T. Kaplan and L. J. Gray, J. Phys . C 9, L483 (1976). 

95. T. Kaplan and L. J. Gray, Phys. Rev. B 15, 3260 (1977). 

96. S. S. A. Razee and R. Prasad (unpublished). 

97. P. R. Halmos, Measure Theory (Van Nostrand, New York, 
1950 ) . 

98. W. Feller, An Introduction to Probability Theory and its 
Applications (Wiley, New York, 1968). 

99. J. W. D. Connoly and A. R. Williams, Phys. Rev. B 27, 
5169 (1983). 

100. For a review on alloy phase stability see D. de Fontaine, 
in Electronic Band Struc ture and its Applications , edited 
by M. Yussouff ( Springer-Verlag , Berlin, 1987), p. 410. 

101. H. Dreysse, A. Berera, L. T. Wille, and D. de Fontaine, 
Phys. Rev. B 39 , 2442 (1989). 

102. P. E. A. Turchi, G. M. Stocks, W. H. Butler, D. M. 
Nicholson, and A. Gonis, Phys. Rev. B 37, 5982 (1988). 

103. For a comprehensive review on disordered systems see J. 
M. Ziman, Models of Disorder (Cambridge University Press, 
Cambridge, 1982). 

104. J. S. Faulkner, Phys. Rev. B 1, 934 (1970). 

105. L. M. Roth, Phys. Rev. B9, 2476 (1975); Phys. Rev. B 
11 , 3769 (1975) . 

106. L. Huisman, D. Nicholson, L. M. Schwartz, and A. 

Phys, Rev. B 24, 1824 (1981). 


Bansil , 



. * * 


172 


107. L. M. Schwartz and H. Ehrenreich, Ann. Phys. 64, 100 

(1971 ) , 

108. B. L. Gybrffy, Phys. Rev. Bl, 3290 (1970). 

109. J. Korringa and R. Mills, Phys. Rev. B 5, 1654 (1972). 
no. L. M. Schwartz, Phys. Rev. B 21 , 522 (1980). 

111. P. Weinberger, R. Dirl, A. M. Boring, A. Gonis, and A, 

J- Freeman, Phys. Rev. B 37, 1383 (1988) and references 

therein . 

112. J. Banhart, P. Weinberger, and J. Voitlander, Phys. Rev. 
B 40, 12079 (1989 ) . 


